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PREFACE 


This  publication  (Interim  Report  69-3)  represents  to  a  large  degree 
the  culmination  over  ten  years  of  research  to  understand  the  integrated 
soil-plant-atmosphere  system.  It  lias  been  a  team  effort  of  technicians, 
professionals  and  graduate  students.  In  the  beginning,  it  was  necessary 
to  study  separately  the  component  part  physical  and  physiological  processes 
under  field  conditions.  This  necessitated  development  of  suitable  measuring 
and  data-handling  techniques.  Once  there  was  sufficient  knowledge  and 
expertise  the  team  was  able  to  study  simultaneously  the  component  parts 
all  together.  Finally,  Douglas  Stewart  developed  a  suitable  computer  model 
to  simulate  a  plant  community  —  both  its  environment  and  its  interaction  — 
then  tested  the  model  against  some  of  the  integrated  field  measurements. 

The  model  is  not  perfect.  It  shouldn't  be,  for  we  do  not  unders'^nd  all 
the  physics  and  physiology. 

Early  in  1970  two  more  Interim  Reports  will  be  published.  One  deals 
with  extensive  testing  of  the  model  with  emphasis  on  plant  water  relations. 
The  other  report  takes  up  in  detail  the  field  instrumentation  used  in 
gathering  the  test  data. 

There  follows  a  listing  of  the  Interim  Reports  published  to  date  by 
the  team. 


E.  R.  Lemon 


i 

f 


LIST  OF  INTERIM  REPORTS 


< 


Interim  Report 


»  62-1 

Mathematical  Study  of  the  First  Stage  of  Drying  of  a  Moist  Soil 
—  W.  Covey 

*  62-2 

Theoretical  Estimate  of  Photosynthesis  under  Field  Conditions  — 

C.  S.  Yocum 

,  62-3 

Solar  Radiation  Balance  and  Photosynthetic  Efficiency  —  C.  S.  Yocum, 

L.  H.  Allen  and  E.  R.  Lemon 

*1 

*62-4 

Radiant  Energy  Exchanges  within  a  Corn  Crop  Canopy  and  Implications 
in  Water  Use  Efficiency  —  L.  H.  Allen,  C.  S.  Yocum  and  E.  R.  Lemon 

62-5 

Analysis  of  Micrometeoro logical  Vertical  Profiles  for  the  Evaluation 
of  Surface  Characteristics  and  Vertical  Fluxes  —  W.  Covey  ' 

„  62-6 

Soil  Moisture  Analyses,  Ellis  Hollow  (Ithaca,  N.Y.),  1960  and  1961 
—  J.  H.  Shinn,  K.  W.  Brown  and  R.  F.  West 

ii 


TABLE  OF  CONTENTS 

Chapter  Page 

INTRODUCTION .  1 

I.  RADIATION  MODELS .  4 

Review . 5 

Theoretical  Considerations  .  5 

Visible  Radiation  .  7 

Near  Infrared  Radiation  .  7 

Thermal  Radiation . . .  l6 

All  Components  . .  18 

Materials  and  Methods  .  .  1 .  19 

Light  Measurements .  20 

Leaf  Measurements .  26 

Results  and  Discussions  .  29 

Conclusions . 45 

II.  A  LEAF  PHOTOSYNTHESIS  MODEL .  46 

Theoretical  Considerations  .  hj 

Materials  and  Methods .  52 

Results  and  Discussions  .  55 

Conclusions .  6l 

III.  THE  BOUNDARY  LAYER  RESISTANCE  OF  PLANT  LEAVES  .  6 2 


Theoretical  Considerations 


65 


Chapter  Page 

Materials  and  Methods  .  68 

Results  and  Discussions  .  TO 

Conclusions . 78 

IV.  SIMULATING  THE  CROP  RESPONSE . 79 

Theoretical  Considerations  ....  .  8l 

Above  the  Vegetation . 81 

Exchange  Within  the  Crop . 87 

Sensible  Heat  and  Mass  Exchange . 89 

Leaf  Temperatures . 97 

The  Soil  Surface  .  . . 99 

Successive  Approximations  .  101 

Materials  and  Methods  .  103 

Results  and  Discussions  .  110 

Conclusions  . . 125 


v 


LIST  OF  TABLES 


Table  Page 

I  Temperature  Coefficients  .  2k 

II  Mean  theoretical  and  energy  balance  values  for  the  total 

flux  of  sensible  and  latent  heat  and  of  CO2  from  the 
unthinned  corn  crop .  122 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1.  Percent  transmission  through  a  hypothetical  array  of  leaves 
as  a  function  of  the  number  of  increments  into  which  the 

array  is  divided .  11 

2.  Leaf  plane  at  an  angle  IL  to  the  horizontal .  12 

3.  Scattered  Radiation  .  17 

k.  The  custom  mounted  visible  radiation  sensor .  22 

5.  The  percentage  error  from  the  cosine  response  for  three  individual 

cells .  23 

6.  A  leaf  segment  (x^y^),  (x2,y2)  between  1  and  3  dm  of 

height .  28 

7.  Response  of  sensors  as  a  function  of  distance  along  the  sampling 

span  in  a  stand  of  com .  30 

8.  Cumulative  leaf  area  index  (CUM.  L.A.I.)  and  leaf  area  density 

for  two  stands  of  com  in  1968  .  31 

9.  Cumulative  leaf  angle  distribution  for  the  stand  o,f  com  of 

L.A.I.  =  2.60  (Study  made  in  1968) . \  32 

10.  Cumulative  leaf  angle  distribution  for  the  stand  of  corn  of 

L.A.I.  =  3.63  (Study  made  in  1968) .  .33 

11.  Cumulative  leaf  area  index  (CUM.  L.A.I.)  for  the  stand  of  corn 

in  1967  .  3h 

12.  Cumulative  leaf  angle  distribution  for  the  stand  of  corn  of 

5.2  (Study  made  in  1967) . .  .  35 

13.  Percent  transmission  as  a  function  of  cumulative  leaf  area 

index  (CUM.  L.A.I.)  DeWit  and  Duncan  Model .  37 

ll+.  Percent  transmission  as  a  function  of  cumulative  leaf  area 

index  (1967  studies) .  38 

15.  Percent  transmission  as  a  function  of  cumulative  leaf  area 
index  at  various  sun  angles.  (1968  study,  L.A.I.  =  3.63  and 
2.60) . 39 


Figure  Page 

1 6.  Percent  transmission  as  a  function  of  cumulative  leaf  area 

index  (1967  study,  L.A.I.  =  3.13  and  2.38) .  1+0 

17.  Percent  transmission  as  a  function  of  cumulative  leaf  area 

index  when  the  sky  was  completely  overcast .  1+2 

18.  Reflectance  of  visible  radiation  from  the  crop  as  a  function 

of  sun  angle  . . 1+3 

19.  The  resistance  network  for  a  leaf:  a  modification  of  the  net¬ 
work  by  Lake  (19^7) .  1+8 

20.  A  simplified  resistance  network  for  a  leaf .  1+8 


21.  Net  photosynthesis  (N)  as  a  function  of  light  flux  density  (i) 

at  a  (#2  concentration  of  260  ppm  (Measurements  only)  ....  56 

22.  Net  photosynthesis  (N)  as  a  function  of  light  flux  density  (i) 

at  a  CO2  concentration  of  260  ppm  (Measurements  and  theory)  .  58 

23.  Net  photosynthesis  (N)  as  a  function  of  CO2  concentration  at 


a  light  flux  density  of  0.230  p  einsteins/cn^/sec  .  59 

2l+.  Boundary  layer  resistance  for  heat  (ra^)  as  a  function  of 

(L,/u)1'2(Pr)  for  a  leaf  in  the  horizontal  position .  71 

25.  Boundary  layer  resistance  for  heat  (ra^)  as  a  function  of 

(L1 /u)1' 2(Pr)  for  a  leaf  in  an  upright  position .  73 

26.  Boundary  layer  resistance  for  heat  (ra^)  as  a  function  of 
(L'/i)1'2  (Pr)  for  a  leaf  in  an  upright  position  with  two 

leaves  upwind .  7I+ 

27.  Boundary  layer  resistance  for  heat  (ra^)  as  a  function  of 

( L*  /u)  (Windtunnel  and  field  measurements) .  76 

28.  The  ratio  of  the  eddy  diffusivity  of  heat  (KH)  and  momentum  (%) 
as  a  function  of  the  instability  parameter  (Z/L)  after 

Panofsky  (1969)  ..........  .  85 

29.  A  represented  of  the  crop  divided  into  n  increments .  91 

30.  Wind  velocity  (u)  as  a  function  of  height  (z) . 110 

31(a).  Eddy  diffusivity  for  sensible  heat  (%)  as  a  function  of 

height  (z)  (Aug.  15,  1200) . Ill 

31(b).  Eddy  diffusivity  for  sensible  heat  (Kh)  as  a  function  of 

height  (z)  (Aug.  18,  1200) .  112 


Figure  r  Page 

31(c).  Eddy  diffusivity  for  sensible  heat  (IC,)  as  a  function  of 

height  (z)  (Aug.  18,  0900) . 7 .  113 

32(a).  Theoretical  profile  of  COp  and  water  vapor  concentration  with 

measurements  (Aug.  15,  12&0) .  116 

32(b).  Theoretical  profiles  of  air  temperature  and  net  radiation  with 

measurements  (Aug.  15,  1200)  .  117 

33(a).  Theoretical  profiles  of  CO2  and  water  vapor  concentration  with 

measurements  (Aug.  18,  1200)  .  118 

33(b).  Theoretical  profiles  of  air  temperature  and  net  radiation  with 

measurements  (Aug.  18,  1200)  .  .' .  119 

3U(a).  Theoretical  profiles  of  CO2  and  water  vapor  concentrations  with 

measurements  (Aug.  18,  0900 ) .  120 


3^(b).  Theoretical  profiles  of  air  temperature  and  net  radiation  with 
measurements  (Aug.  18,  0900 ) . . 


121 


'-jv  '  ■  ■. 


INTRODUCTION 


In  general,  the  age  old  question  of  what  constitutes  crop  yield 
lies,  as  yet,  unanswered.  Hie  plant  when  it  germinates  produces  a  cer¬ 
tain  amount  of  leaf  area  from  sources  of  carbohydrate  stored  in  the 
seed.  It  uses  this  leaf  area  to  intercept  radiation  which  it,  in  turn, 
uses  to  convert  CO2  to  carbohydrate.  The  carbohydrate  provides  energy 
and  building  material  for  the  creation  of  roots,  stem  and  more  leaf  area. 
While  simply  stated  here,  the  processes  of  germination,  vegetative  growth, 
flowering  and  fruiting  are  extremely  complex  and  include  many  inter¬ 
actions  with  a  constantly  changing  environment.  Because  of  these  com¬ 
plexities,  the  traditional  methods  of  growth  analysis  described  by 
Watson  (19^7,  1952),  Williams  (19^6)  and  others  along  with  correspond¬ 
ing  statistical  techniques  such  as  the  studies  of  Black  (1955)  and 
Glenday  (1956)  have  not  provided  an  understanding  of  the  mechanisms  by 
which  a  plant  is  influenced  by  and  in  turn  responds  to  the  environment. 

An  alternative  to  statistical  correlation  is  the  method  of  mathe¬ 
matically  modelling  a  plant-environment  system.  That  is,  the  various 
growth  processes  are  described  in  mathematical  terms  and  integrated 
with  time.  Davidson  and  Philip  (1956)  presented  one  of  the  first  simu¬ 
lation  models  where  they  considered  only  the  light  environment  and  its 
effect  on  the  accumulation  of  dry  matter.  Since  then  computers  with 
their  ability  to  handle  detail  have  been  used  in  not  only  complex 
radiation  models  (de  Wit  ( 19 65),  Duncan,  et_  al .  (1967))  but  also  in 
growth  simulation  models  (Brower  and  de  Wit,  1968).  Hie  latter  calculates 
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the  daily  accumulation  of  net  photosynthesis  and  then  models  the  trans¬ 
location  processes  vhere  carbohydrates  are  distributed  to  the  stem  and 
root  tissues. 

Although  models  are  available  which  describe  the  light  regime  in  a 
plant  system  in  detail,  less  attention  has  been  applied  to  the  problem 
of  COg  exchange  and  the  dissipation  of  radiation  loads  on  individual 
leaves  through  sensible  and  latent  heat  exchange.  In  particular,  under¬ 
standing  of  the  effects  of  the  architecture  of  the  plant  system  on  these 
exchange  processes  is  lacking. 

There  are  two  types  of  field  studies  which  are  essential  to  the 
development  of  a  complete  model  of  the  plant-environment  system.  The 
first  is  the  investigation  of  net  assimilation  of  C02  of  individual 
leaves  using  leaf  chambers.  (For  example,  Gaastra  (1959),  Hesketh  and 
Musgrave  (1962)).  In  particular,  net  photosynthesis  of  individual 

leaves  must  be  known  as  functions  of  ligjht  flux  density,  CCU  concentration 

/ 

at  the  leaf  surface  and  leaf  temperature  (effects  of  water  stress  will 
not  be  considered  directly  in  this  treatment).  These  functions  are  used 
directly  in  the  simulation  model  described  in  this  study. 

The  second  type  of  studies  necessary  for  simulation  are  the  in  situ 
investigations  made  by  agrometeorologists.  (Lemon  (i960),  Monteith  and 
Szeicz  (i960)).  In  these  studies  profiles  of  air  temperature,  wind 
velocity,  COp  and  water  vapor  concentrations  and  net  radiation  are 
measured  above  and  within  a  plant  community  disturbing  it  as  little  as 
possible.  Aerodynamic  and  energy  balance  methods  using  these  profiles  pro 
vide  estimates  of  net  fluxes  of  C02  and  sensible  and  latent  heat  (Lemon 
1965).  Thus  these  fluxes  can  be  compared  to  simulated  values  since  the 


calculations  are  independent  of  the  above  methods.  Also  profiles  can  be 
generated  by  the  simulation  model  and  compared  with  measurements.  Finally 
the  above  measurements  made  2-3  meters  above  the  crop  and  at  the  soil 
surface  provided  needed  boundary  values  for  the  mathematical  solutions. 

To  construct  this  simulation  mo^el,  a  rather  wide  range  of  informa¬ 
tion  must  be  assembled.  First  of  all,  radiation  loads  on  individual 
leaves  in  a  plant  community  must  be  known.  A  procedure  will  be  outlined 
to  divide  the  leaf  area  of  a  crop  into  classes  of  radiation  loads.  This 
is  done  for  visible  radiation  in  conjunction  with  the  calculation  of  net 
photosynthesis  as  well  as  total  radiation  loads  used  in  estimating  sensi¬ 
ble  and  latent  heat.  Because  of  the  immense  amount  of  calculations  in¬ 
volved,  a  computer  is  used.  Thus  in  subsequent  theoretical  discussions, 
programming  techniques  are  introduced  when  they  were  considered  essential 
to  the  solution  of  the  equations  involved. 

Once  the  distribution  of  radiation  in  a  crop  system  is  known,  the 
photosynthetic  response  of  individual  leaves  to  radiation  can  be  inte¬ 
grated  over  the  complete  crop  cover.  However,  CO2  supply  and  leaf 
temperatures  also  influence  these  responses.  Thus  the  exchange  processes 
between  the  individual  leaves  and  the  air  stream  both  above  and  within 
the  crop  have  to  be  introduced  in  this  integration  scheme.  Since  leaf 
responses  and  properties  of  the  air  stream  are  interdependent,  the  equa¬ 
tions  describing  the  exchange  processes  of  the  air  flow  and  the  equations 
determining  the  leaf  responses  are  solved  simultaneously  using  success¬ 
ive  approximation  techniques.  Here  again  a  computer  was  used  extensively 
in  these  solutions .  The  development  of  the  various  components  of  the 
integration  scheme  along  with  experinental  testing  are  considered  in 


subsequent  chapters. 


Chapter  I 
RADIATION  MODELS 

As  stated  in  the  introduction  the  radiation  regime  in  the  crop 
is  divided  into  visible,  near  infrared  and  thermal  regions  of  the  spectrum. 
Visible  radiation  supplies  the  excitation  energy  to  the  chlorophyll 
molecule  for  photosynthesis.  It  also  determines  to  a  large  extent  the 
width  of  stomatal  apertures  when  the  plant  is  not  under  water  stress. 

Thus  the  emphasis  here  both  theoretically  and  experimentally  will  be 
on  visible  radiation.  Theoretical  treatments  of  the  near  infrared  and 
thermal  components  are  also  given. 
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Review 

Visible  radiation  in  crop  systems  has  been  extensively  studied  mainly 
because  of  its  direct  effect  on  photosynthesis.  Anderson's  review  (?96U) 
provides  some  indication  of  the  large  interest  in  this  area,  therefore, 

I  will  not  attempt  a  complete  review  of  this  subject  but  will  cover  only 
the  studies  which  are  pertinent  to  the  theory  developed  later  in  this 
discussion. 

The  most  widely  used  equation  describing  light  penetration  is  the 
following  analogy  of  Beer's  Law; 

I  =  IQ  e*^  (1) 

F  is  the  cumulative  leaf  area  per  unit  soil  area  measured  from  the  top 
of  the  crop  downward; 

K  is  an  extinction  coefficient; 

I  is  the  flux  density  of  visible  radiation  at  F; 

and  IQ  is  the  corresponding  flux  density  at  the  surface  of  the  crop. 

Monsi  and  Saeki  (1953),  Takeda  and  Kumera  (1957)  and  others  found  that 
this  equation  adequately  described  light  penetration  in  many  types  of 
plant  communities.  Davidson  and  Philip  (1958)  used  equation  (l)  (sett¬ 
ing  K=l)  in  their  growth  simulation  model. 

Since  then,  Isobe  (1962),  Anderson  (1966)  and  Chartier  (1966a)  among 
others  have  derived  equation  (l)  by  assuming  a  random  distribution  of 
leaf  area.  K  is  a  property  of  the  mean  leaf  angle  and  therefore  may  or 
may  not  be  constant  with  F. 

Other  equations  for  light  penetration  have  been  developed  for  diff¬ 
erent  distributions  of  leaf  area.  For  example,  Monteith  ( 19 65)  used  a 
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positive  binomial  distribution;  Niilisk,  et  al .  (1969 )  used  a  negative 
binomial  and  a  Markov  distribution. 

f  A  part  of  the  radiation  fall ing  on  plant  leaves  is  transmitted  or 
reflected.  Kasanaga  and  Monsi  (195*0  extended  the  exponential  law  (Equa¬ 
tion  l)  to  include  the  effects  of  leaf  transmission.  A  different  approach 
considers  the  crop  as  a  homogeneous  system  characterized  by  a  transmission 
coefficient  (closely  related  to  K)  and  a  reflection  coefficient.  The 
Kubelka-Monk  theory  presented  in  detail  by  Allen  and  Richardson  (1968) 
does  this.  This  theory  was  compared  with  measurements  by  Allen  and 
Brown  (1965).  In  their  crop  system,  the  transmission  coefficient  did 
not  remain  constant  with  F. 

There  is  a  considerable  advantage  to  study  the  light  regime  in  a 
canopy  in  terms  of  the  characteristics  of  individual  leaves  rather  than 
the  crop  as  a  whole.  This  allows  one  to  take  into  account  differences 
in  leaf  properties  with  depth  and  eventually,  differences  in  leaf  area 

i 

distributions  with  depth.  There  have  been  at  least  two  attempts 
(de  Wit  (1965)  and  Duncan  et  al.  (196?))  which  use  leaf  angle  distri¬ 
butions  and  leaf  transmission  and  reflection  coefficients  in  describing 
the  light  regime  in  plant  communities.  These  studies  are  particularly 
important  since  the  light  irradiance  levels  crossing  leaf  plants  in  the 
crop  are  calculated.  The  following  theoretical  considerations  rely 


heavily  on  these  studies. 
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Theoretical  Considerations 

The  three  regions  of  the  spectral  band  of  radiation  considered 
here,  have  quite  different  leaf  transmission  and  reflection  coefficients. 
For  visible,  they  are  of  the  order  of  0.1;  for  near- infrared  they  vary  from 
0.4  to  0.5  and  for  long  wave  radiation  they  are  approximately  zero.  These 
coefficients  vary  with  wavelength  to  some  extent  within  these  regions 
but  this  variation  is  small  when  compared  with  the  differences  between 
regions.  The  theory  dealing  with  each  region  will  be  considered  in  turn. 

Visible  Radiation 

De  Wit  (1965)  introduced  the  following  equations  relating  leaf 
angle  (IL),  sun  angle  (IS),  the  angle  between  the  sun  and  the  leaf  (LS) 
and  the  difference  between  the  azimuths  of  the  sun  and  leaf  planes  (DA). 


Sin(LS)  =  A  + 

B  sin (DA) 

(2) 

A  =  sin (IS)  ♦ 

cos(IL) 

(3) 

B  =  cos (IS)  • 

sin(lL) 

(4) 

Sin(LS)  is  also  the  projection  of  the  leaf  in  the  direction  of  the  sun's 
rays.  Sin(LS)  divided  by  sin(lS)  is  the  projection  of  a  unit  leaf  area 
onto  the  horizontal  and  is  directly  related  to  the  amount  of  light 
intercepted  by  the  leaf.  To  obtain  a  mean  projection  (PL),  one  must 
integrate  equation  (2)  from  -90°  to  +90°. 
w/2 

/  (A  +  B  sin (DA))  d(DA) 

— tt/2 


That  is 


PL 


/  d(DA) 
—tt/2 


(5) 


If  IS  >  IL 


PL  =  A 


Wipf 


_ : — iL. - 
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If  IS  <  IL  PL  =  2/tt [B  cos ( DAO )  -  A(DAO)  ]  (7) 

where  DAO  =  arcsin  (-A/B) 

Defining 

0o  =  ( tt/2  -  DAO) 
equation  (7)  becomes 

PL  =  2/n[sin0o  cos(lS)  sin(lL)]  +  (l  -  0o/9O°)  cos(IL)  sin(IS)  (8) 

Thus  PL  is  identical  to  F/F'  defined  in  appendix  A  of  Duncan,  et  al. 
(1967)*  For  a  given  increment  of  leaf  area,  the  effects  of  the  various 
leaf  angles  have  to  be  included,.  Thus  the  mean  projection  (PL)  of  the 
leaf  area  increment  becomes 

_  9 

PL  =  E  PL,  *  FR 

J-l  J  J 

where  FR^  is  the  relative  amount  of  leaf  area  in  the  leaf  angle  class  j . 
Nine  leaf  angle  classes  are  used  in  this  study. 

In  studying  light  penetration,  a  crop  can  be  divided  into  a  number 
of  finite  increments  with  a  leaf  area  of  S  in  each  increment.  The  per¬ 
cent  transmission  (PT)  of  light  to  but  not  including  the  n+1  layer  is 
given  by  de  Wit  (1965)  as; 

PT  =  100  (1  -  S  »  (PL)  )n  (9) 

sin( IS) 

Duncan,  et  al .  (1967)  derived  a  different  expression  using  equations  by 
Wilson  (i960).  ,The  mean  probability  M  of  a  sun’s  ray  being  intercepted 
by  a  leaf  was  defined  as 

M  -  s  : 

"  sin (IS)  ' 
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If  leaves  are  distributed  randomly,  percent  transmission  can  be 
expressed  as: 

PT  =  100(1  -  [2*2.”“  +  +  Mk"M  +....])  (10) 

1!  2!  3! 

The  terms  inside  the  square  brackets  sure  probabilities  of  1,  2,  3  etc. 
leaves  being  in  the  path  of  a  sun's  ray.  Equation  (10 )  reduces  to 

PT  =  100(e'M)n  (11) 

Note  that  equation  (8)  is  a  function  of  S  while  (ll)  is  not.  Figure  (l) 
shows  percent  transmission  of  a  hypothetical  array  of  leaves  as  a  func¬ 
tion  of  the  number  of  increments  into  which  the  leaf  area  was  divided. 

De  Wit  suggested  that  S  be  used  as  an  empirical  parameter  to  fit  measured 
light  profiles  with  theoretical  curves.  However  for  S  less  than  0.4, 
small  changes  occur  in  PT  using  de  Wit's  equation.  Note  that  as  n  becomes 
larger  equation  (8)  and  (ll)  give  approximately  the  same  answer.  Also 
note  that  in  equation  (8)  and  (ll),  PL  appears  sis  a  constant.  In  the 
actual  computer  program,  PL  varies  with  the  leaf  angle  distribution 
which  may  vary  with  F. 

Thus  far,  only  penetration  of  direct  radiation  from  the  sun  has-been 
considered.  Sky  radiation  can  be  treated  in  the  same  way  if  we  assume 
the  sky  is  a  uniformly  bright  hemisphere,  (de  Wit,  1965).  The  sky 
is  divided  into  9  angle  classes  from  the  horizon  to  the  zenith.  The 
proportion  of  light  (PD^)  falling  on  the  horizontal  from  a  sky  segment 
from  IS^  to  ISg  is  expressed  as 


2 

/  2ira  cos  (IS)  *  sin(lS)  d(lS) 


tt/2 

/  2ira  cos  (IS)  •  sin(lS)  d(is) 


0 

=  sin2(lS2)  -  sin2(lSq) 


(12) 


NQ  of  INCREMENTS 


Figure  1 


Per  cent  transmission  through  a  hypothetical  array  of  leaves 
as  a  function  of  the  number  of  increments  into  which  the 
array  is  divided.  L.A.I.''-=,  4.0  (is)  =  60°,  IL  =  45°.  The 
straight  line  represents  Equation  (ll);  the  asymptotic  line 
represents  Equation  (9). 


I  "  *  «■  »  * 
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'a'  is  an  imaginary  radius  of  the  sky  hemisphere.  PD^  times  the  total 
amount  of  sky  radiation  is  the  amount  of  sky  radiation  arriving  at  the 
angle  IS^  where  in  this  case 


is.  =  (n^  +  is2)/2. 


Special  attention  is  given  to  the  proportion  of  sky  light  which  pene¬ 
trates  each  increment  leaf  area  of  the  crop.  These  values  are  stored  in 
an  array  TIk  where  i  refers  to  the  increment  of  leaf  area.  These  values 
are  used  repeatedly  in  the  program. 

The  flux  density  of  light  striking  a  given  leaf  in  the  canopy  is 

sin(LS)  •  PI 
sin (IS) 

where  (Dl)  is  the  flux  density  of  direct  radiation  on  a  horizontal  sur¬ 
face.  The  amount  of  light  transmitted  and  reflected  are  respectively 
(sin(LS)  •  (DI)  •  A)/sin(lS)  and  (sin(LS)  •  (Dl)  •  y)/sin(lS),  where  A 
and  y  are  the  leaf  transmissivity  and  reflectivity  respectively. 

Two  simplifying  assumptions  are:  (a)  that  light  is  scattered  by 
the  leaf  diffusely  and  (b)  that  only  two  directions  (upward,  down)  of 
scattered  radiation  need  be  considered.  Figure  (2)  shows  an  element  of 
leaf  dA  at  angle  (IL)  to  the  horizontal.  The  amount  of  light  reflected 
(W)  from  dA  is  expressed  as 

U  =  Sln(IJ3>  •  DI  •  I>  (11.) 

sin(lS) 


W  =  /  27tW0  i^sinB  cos0  d6 

0 


where  0  is  the  angle  of  elevation  from  the  leaf  plane,  W0  is  the  intensity 
of  radiation  emitted  perpendicular  to  the  plane  and  r  is  the  radius  of 
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an  imaginary  hemisphere.  In  figure  (2),  if  ML  =  z ,  KL  =  y,  K0=  x  and 
OM  =  r  then  sin0^  =  z/r,  tan6  =  y/x  and  tan(lL)  =  z/y.  8  is  the  angle 
about  the  normal  to  the  leaf  plane  at  0;  0^  is  the  angle  of  elevation 
bounded  by  the  plane  AMEC. 

Also  x2  +  y2  +  z2  =  r2  -  (l6) 

Combining  the  above  equations,  one  obtains  .  . 


=  arc  sin 


B1  sin20 
1  +  B'sin20 


where  B'  =  tan^IL 

The  amount  of  radiation  V  scattered  from  dA  through  the  portion  of  the 
hemisphere  between  AMEC  and  the  plane  of  the  leaf  is  expressed  as 


W '  =  /  /  W0r  sin0  cos0  d0  d<j> 

0  0 


_  r2  u  f*  B'  sin2  <j> 


-"o/ 

2  0 


1  +  B’sin  $ 


The  fraction  (e)  of  W  to  W  is  given  as 

E  =  W'/W  =  l/(2ir)  /  __B'sin2<t> 

0  1  +  B'sin2<)> 


E  =  1/2  - 


+  2B' 


The  fraction  (E)  is  both  the  proportion  of  radiation  reflected  downward 
and  the  proportion  transmitted  upwards,  (l  -  E)  is  the  converse. 
Therefore  the  total  amount  of  light  scattered  upward  by  the  leaf  (Set) 
and  downward  (Sc+)  are  respectively; 


Set  =  sin(LS)  •  PI 
sin(lS) 


(A(l  -  E)  +  iiE) 


and  Sc+  = 
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(21) 


sin(LS)  •  DI 
sin(lS) 


(AE  +  u(l  -  E)) 


In  the  program,  PL^  replaces  sin(LS)  where  j  represents  a  leaf 
angle  class.  PLj  is  the  light  flux  density  at  the  leaf  surface  of 
angle  (j)  averaged  over  various  positions  around  the  azimuth.  To  obtain 
the  total  amount  of  scattered  radiation.  Set  and  Sc+  are  summed  over  the 
leaf  angle  classes  for  both  the  direct  radiation  and  the  nine  classes 
of  sky  radiation.  A  final  adjustment  is  made  using  equation  (ll). 
Percent  transmission  is  a  measure  of  the  actual  area  of  leaves  reached 
by  direct  and  sky  radiation. 

At  this  stage,  direct  and  sky  radiation  penetrate  into  a  crop 
creating  sources  of  scattered  radiation.  Radiation  is  assumed  to  be 
scattered  equally  in  all  directions.  Thus,  this  scattered  radiation 
appears  as  a  uniformly  bright  infinite  plane  to  other  leaves.  Assuming 
this,  one  can  use  the  procedures  already  worked  out  for  sky  radiation  to 
calculate  the  penetration  both  up  and  down  of  radiation  scattered  by  the 
leaves  in  the  crop.  These  amounts  of  scattered  radiation  which  pene¬ 
trate  to  other  levels  is  simply  added  to  the  previous  values  of  Set  and 
Set.  Also  the  amounts  of  sky  radiation  penetrating  directly  into  the 
canopy  are  added  to  Set.  Thus  one  can  now  define  Sc^t  and  Sc^t  as  the 
total  amounts  of  diffuse  radiation  in  the  crop  moving  up  and  down 
respectively  at  level  i. 

The  soil  surface  in  this  study  was  treated  as  if  it  were  a  layer 
of  leaves  with  A  =  0.0  and  u  =  0.13. 

Of  primary  interest  in  a  photosynthesis  model  is  the  flux  densities 

(Lt.  )  of  visible  radiation  on  individual  leaves  in  the  canopy.  The 
ij 
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equation  defining  the  radiation  load  on  a  leaf  is  given  as 

Lt.  =  I A  +  B  sinDAl  •  -  +  Sc.  +  +  Sc. + 

13  siallS)  1+1 

for  leaves  in  direct  sunlight  and 


(22) 


Ltij  =  S<W  +  Sci  + 


(23) 


for  leaves  shaded  from  the  sun.  Leaves  in  direct  sunlight  were  consid¬ 
ered  at  18  positions  of  (DA)  between  -ir/2  to  +7r/2.  Also  the  relative 
amount  of  leaf  area  A^  in  direct  sunlight  for  a  given  increment  i  was 
found  from 

Aoi  =  (PTj  -  PTi+1)  •  -!la(l°>  (24) 

(PLi  *  100) 

PT^  is  the  percent  transmission  of  direct  radiation  to  the  i^  increment 
(Duncan  et_  al.  1967).  Equations  22-2U  were  used  to  determine  total  visible 
radiation  loads  at  all  leaf  positions  in  the  crop.  These  loads  were 
expressed  in  classes  as  0-5^,  5-10^,  ±0-Vy%  of  Dl/sin(lS).  Leaf  areas 
with  radiation  loads  in  the  same  range  were  summed.  Thus  a  frequency 
distribution  of  leaf  area  in  these  classes  for  each  increment  was  the  only 
information  stored  from  the  above  calculations. 

Another  advantage  to  this  procedure  is  the  ability  to  separate  leaf 
areas  with  direct  radiation  loads  falling  on  the  dorsal  as  opposed  to 
the  ventril  sides  of  leaves.  This  happens  whenever  (A  +  B  sinDA)  is 
negative.  Therefore,  a  separate  frequency  distribution  of  these  leaf 
areas  could  be  constructed  although  in  this  study  this  information  was 
not  used.  Both  sides  of  the  leaf  were  assumed  to  have  equal  responses 


to  light. 
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Near  Infrared  Radiation 

The  near  infrared  radiation  can  be  treated  in  the  same  way  as  the 
visible  radiation  with  one  important  difference.  Multiple  scattering 
of  visible  radiation  was  considered  negligible  because  of  the  relatively 
low  leaf  reflection  and  transmission  coefficients.  Niilisk  et_  al .  (1969) 
made  this  same  assumption.  However  the  relatively  large  values  of  these 
coefficients  in  the  near  infrared  make  a  consideration  of  multiple 
scattering  necessary. 

Radiation  from  the  sun  and  sky  penetrates  to  a  given  leaf  area 
increment  and  creates  first  order  sources  of  scattered  radiation.  Some 
of  this  radiation  reaches  the-  soil  and  other  leaves  and  is  reflected 
back  to  the  same  increment  to  create  second  order  sources.  This  is  illus¬ 
trated  in  Figure  (3).  Let  S^t  be  the  amount  of  light  scattered  upward 
and  Si  the  amount  scattered  downward  at  level  i  from  the  first  order 
sources.  The  amount  (R^l)  of  reflected  back  is  expressed  by 


V  =  Si+ 


p[l  +  (TR  !  +  SX)2  + 


(TR,_,  +  SX)2(TRj_2  +  SX)2  + 


J-l 


Similarly 

EA  *  S^  •  S  •  y [1  +  (TRj  +  SX)2  +  ( TRj  +  SX)2(TRJ+1  +  SX)2  +  . 


(25) 

.  ] 

(26) 


In  this  case  j  represents  the  increment  of  leaf  area  directly  below  level 
i.  The  assumptions  involved  in  equations  (25)  and  (26)  are  that  the 
scattered  radiation  penetrates  a  given  leaf  area  increment  in  he  same 
relative  amounts  as  sky  radiation.  Also,  the  total  amount  of  radiation 
f-eflected  and  transmitted  by  a  leaf  area  increment  is  directly  propor¬ 
tional  to  SX  and  Sp  respectively  where  S  is  the  leaf  area  per  increment. 

The  secondary  sources  then  scatter  radiation  upwards  (SS^t)  and 
downwards  (SS^+)  expressed  respectively  as 


» 


•s 


(2?) 
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SSi+  =  R^  •  Sy  +  Ri+1t  •  SA 

SS^  =  R^  •  SA  +  Ri+1+  •  Sy  (28) 

Let  3+  =  SS^t/S^t  and  6+  =  SS^+/S^+.  Then  the  total  amounts  of  radiation 
scattered  upward  and  downward  from  increment  i  are  expressed  respectively  as 


Scit  =  Si+(l  +  6+  +  6+2  +  6+3  + . )  (29) 

and  Sc±i  =  Si+(l  +  8+  +  6+?  +  6  +  3  + . )  (30) 


Note  that  the  soil  is  considered  again  as  a  layer  of  leaves  with  zero 
transmittance  and  a  reflectance  of  0.1  for  near  infrared  radiation. 

These  total  sources  of  scattered  near  infrared  radiation  can  now  he 
treated  in  the  same  manner  as  the  first  order  sources  of  visible  radiation. 

Thermal  Radiation 

Thermal  radiation  is  considered  here  for  completeness.  Since 
thermal  radiation  varies  with  leaf  temperatures,  it  is  closely  related 
to  the  energy  budget  equations  which  will  be  considered  in  subsequent 
chapters.  At  this  time,  I  will  assume  all  the  leaf  temperatures  are  known. 

A  leaf  in  the  community  'sees'  'in  an  upward  direction  a  fraction 
of  sky  and  other  leaves.  'Looking'  downward  it  'sees'  a  fraction  of  the 
soil  surface  and  other  leaves,  (idso  (1968)).  The  thermal  radiation 
moving  downward  in  the  crop  TH^  is  expressed  as; 

TH.  =  (F  0  S  Rt)  •  SK  +  (1  -  F  0  S  R)  •  TL.  ( 3l) 

where  F  0  S  R  is  the  fraction  of  thermal  sky  radiation  which  reaches 
the  level  i.  Since  thermal  sky  radiation  (SK)  is  regarded  as  coming  from 


an  infinite  hemisphere  of  uniform  temperature,  this  fraction  is  identical 


to  the  fraction  of  visible  sky  radiation  penetrating  to  level  i  which 
was  calculated  above.  TL^  is  an  average  leaf  temperature  above  the  level 
i.  A  procedure  has  been  outlined  to  divide  leaf  area  into  classes 
based  on  visible  radiation  loads.  A  similar  procedure  can  be  used  to 
divide  the  leaf  area  of  each  increment  into  classes  of  total  solar 
radiation  loads .  In  Chapter  IV  a  procedure  is  outlined  to  calculate 
the  leaf  temperature  of  each  solar*  radiation  class  in  each  increment. 

Thus  TL  is  weighted  according  to  the  leaf  area  in  each  radiation  class 
as  well  as  the  proximity  of  the  leaf  area  to  the  level  i.  Thus  leaf 
temperatures  in  increment  1  would  be  weighted  according  to  the  fraction 
of  diffuse  radiation  which  would  penetrate  from  increment  1  to  increment 
i.  The  stored  values  of  TR  are  used  for  this  calculation. 

A  similar  procedure  can  be  used  to  calculate  TH+  or  upward  directed 
thermal  radiation.  The  amount  of  thermal  radiation  (based  on  a  soil 
temperature)  coming  from  the  soil  replaces  the  sky  radiation.  A  procedure 
to  calculate  soil  temperature  is  given  in  Chapter  IV. 

All  Components 

The  net  radiation  can  now  be  calculated  for  a  given  level  since  all 
components  of  radiation  moving  upward  or  down  can  be  calculated.  Thus; 

Rn.  =  TH.+  -  TH-jt  +  VR+  -  VR+  +  IR+  -  IR+ 

Tli  i  1 

where  VRt  and  VR-t-  are  the  mean  values  of  visible  radiation  moving  upward 
and  downward  respectively.  Similarly  IRt  and  IR  -  are  the  near  infrared 
components.  Thermal  radiation  and  leaf  temperatures  are  discussed  in  more 


detail  in  Chapter  IV.  Also,  comparisons  between  measured  and  theoretical 
profiles  of  net  radiation  are  given  in  Chapter  IV. 
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Materials  and  Methods 

A  verification  of  the  visible  radiation  model  was  attempted  by 
measuring  mean  radiation  flux  densities  at  various  levels  in  stands  of 
com.  These  mean  values  (MV^ )  are  expressed  by 

MV.  =  PTf  (DI  '/(sin(lS)  x  100 )  +  Sc-^  (32) 

DI/sin(lS)  was  measured  directly.  Also  measured  was  Sc^t,  which  is  the 
reflected  radiation  at  the  crop  surface. 

Measurements  were  made  over  a  two-year  period  (1967,  1968)  at  the 
Ellis  Hollow,  Microclimate  Investigations,  Experimental  Site  near  Ithaca, 
New  York.  A  200  acre  field  of  com  was  planted  mechanically  with  a  check 
row  planter  in  an  approximation  of  a  hexagonal  pattern. 

Plants  were  in  38  cm  north-south  rows  with  approximately  U6  cm 
between  plant  in  the  rows.  A  density  of  about  6  plants /m2  was  initially 
obtained.  In  the  1968  study,  one-half  of  the  20  acre  field  was  thinned 
at  emergence  to  U.3  plants/m2.  A  leaf  area  index  of  5.2  was  attained 
in  1967*  In  1968  the  LAI's  varied  depending  on  thinning  and  stage  of 
growth  from  2.56  to  3.63.  In  1968  light  was  measured  at  two  different 
stages  of  growth,  one  stage  Just  before  anthesis  and  the  other  two  weeks 
after  anthesis. 

Li  (jit  Measurements 

The  method  of  measuring  visible  radiation  was  essentially  that  of 
Federer  and  Tanner  (1966b).  A  selenium  cell  was  mounted  under  a  cosine 
correcting  head  and  an  85  c  Wratten  filter.  This  combination  gives  the 
instrument  an  approximation  to  the  photon  response.  That  is,  the  instru¬ 
ment  is  more  sensitive  in  the  red  region  of  the  visible  band  compared  to 


the  blue  and  compensates  for  the  higher  energy  of  the  blue  photon.  Hie 
reasoning  used  here  is  that  if  each  photon  regardless  of  its  wavelength 
contributes  the  same  amount  of  excitation  energy  to  the  chlorophyll 
molecule,  then  a  set  number  of  photons  should  be  the  basic  unit  of 
radiation  rather  than  an  energy  unit.  A  recent  action  spectra  by 
Bulley,  et  al.  (1969)  working  with  radish  leaves  showed  a  variation  of 
±15%  for  constant  incident  quanta. 

In  this  study  each  photocell  was  mounted  in  a  threaded  plexiglass 
holder  which  could  be  screwed  into  a  cosine  correcting  head  as  shown  in 
Figure  (4).  The  head,  machined  from  black  plexiglass  (Kaufman  Glass  Co., 
#2025  Black,  Wilmington,  Delaware),  was  similar  in  shape  to  the  modified 
Scripps  Head  used  by  Kerr,  et  al.  (1967). 

Cosine  correction  depended  on  the  exposed  length  of  the, diffusing 
filter  and  the  distance  from  the  diffusing  filter  to  the  top  inside  edge 
of  the  black  plexiglass .  These  distances ,  were  found  by  trial  and  error 
to  be  0.130  cm  and  0.470  cm  respectively.  The  cosine  response  was  measured 
by  shining  a  collimated  beam  of  light  directly  at  the  diffusing  filter. 

The  instrument  was  mounted  vertically  on  a  rotating  stand.  The  response 
was  measured  as  the  instrument  was  turned  through  +90°  to  -90°.  Per  cent 
deviation  from  the  cosine  law  is  shown  in  Figure  (5)  for  three  different 
instruments . 

The  diffuse  filter  was  made  from  white  acrylic  plastic  (W  =  2447, 
Kaufman  Glass  Co.,  Wilmington,  Delaware).  It  was  needed  not  only  for 
cosine  correction  but  also  to  reduce  the  liggit  intensity.  To  insure  a 
linear  response  to  light  intensity  to  improve  stability  of  the  cell,  a 
.883  cm  length  of  filter  was  used  to  reduce  the  incident  flux  density  to 
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Figure  5.  The  percentage  error  from  the  cosine  response  for  3  individual 


about  8 %  of  the  original.  The  Wratten  85C  filter  (Eastman  Kodak  Co., 
Rochester,  New  York)  was  mounted  underneath  the  diffuse  filter  and  also 


helped  to  reduce  the  light  flux  density  on  the  cell. 

A  precision  carbon  film  resistor  (80.2  dims)  was  mounted  inside  a 
microphone  connector  (amphenol  80M  C  2M).  The  cell  output  was  recorded, 
as  a  change  of  potential  across  the  resistor,  on  a  digital  recorder  (A.  D. 
Data  Systems,  Rochester,  New  York)  accurate  to  ±20  microvolts.  The 
response  of  the  instrument  in  full  sunlight  was  1-2  millivolts. 

There  was  an  advantage  to  using  a  carbon  film  resistor  and  placing 
it  near  the  selenium  cell.  Table  (l)  shows  temperature  coefficients  with 
the  resistor  mounted  both  in  the  instrument  and  outside  of  it. 


TABLE  1 


Temperature  Coefficients 


Unit  with  mounted  resistor 


.00208 


Unit  with  resistor  at  room  temperature 


.00U0T 


The  coefficients  were  measured  by  inserting  a  thermocouple  junction  into 
the  cell  holder  so  that  it  touched  the  edge  of  the  ceT"  without  block¬ 
ing  any  radiation.  The  instrument  was  then  placed  under  a  heat  lamp  of 
constant  intensity  and  its  millivolt  signal  recorded  simultaneously  with 
the  temperature  of  the  cell.  The  coefficients  were  averaged  over  a  change 
of  cell  temperature  from  15°  to  60°C.  The  resistance  of  a  carbon  film 


resistor  decreases  as  its  temperature  increases.  This  tends  to  offset 
„  ..  _ _  •  ...  -  •  *  -  -  *  "•  '  --  • 

the  increase  with  temperature  in  output  from  the  cells. 

In  order  to  obtain  an  average  light  irradiance  at  a  given  level  in 
the  crop  an  adequate  number  of  measurements  have  to  be  made  to  be 
representative  of  that  level.  During  the  course  of  this  study,  two 
sampling  systems  were  used.  Die  first  system,  used  in  1967,  consisted  of 
aluminum  channels  (0.025*+  x  0.025*+  x  5-2  meters)  with  aluminum  sensor 
holders  machined  to  slide  easily  in  the  channels.  Each  light  sensor  was 
attached  to  an  aluminum  holder.  This  holder  was  pulled  from  one  end  of 
the  channel  to  the  other  by  a  pulley  system  attached  to  a  reversible 
gear  motor  (W.  W.  Grainger,  Inc.,  Syracuse,  New  York).  A  tripping 
mechanism  reversed  the  motor  when  a  holder  reached  the  end  of  its  traverse. 

In  1968,  the  aluminum  channels  were  replaced  by  wire  cable.  At  a 
particular  level,  a  selenium  cell  was  placed  in  a  plexiglass  holder 
through  which  two  small  outrigger  cables  were  threaded.  These  outrigger 
cables  were  stretched  between  two  welded  conduit  masts  so  that  they  were 
about  15  cm  apart  in  the  horizontal  plane.  The  masts  were  held  in  place 
by  guy  wires  attached  to  anchor  plates.  Very  high  tensions  were  applied 
to  the  outrigger  cables  and  a  span  of  9  meters  was  obtained  without 
serious  drooping.  The  same  pulley  system,  reversible  gear  motor  and 
tripping  mechanisms  were  used  to  drive  the  cell  back  and  forth  along  the 
span.  The  cells  moved  at  approximately  1.2  meters  per  minute.  Readings 
were  made  every  second.  Five  levels  in  each  stand  of  com  were  sampled. 

Four  larger  versions  of  the  light  sensor  described  above  were  made 
using  a  larger  cell  (Type  A-5,  International  Rectifier  Corp.).  They 
were  used  to  measure  incoming  total  visible,  incoming  diffuse  visible 


(sky)  and  reflected  visible  from  crop  surfaces.  Diffuse  visible  was  measured 
by  placing  a  shade  ring  over  a  sensor  blocking  direct  radiation  from  the 
sun. 

All  the  visible  radiation  sensors  were  calibrated  against  two  Epply 
pyranometers  on  clear  days  using  the  sun  and  sky  as  sources  of  radiation. 

A  standard  single  domed  Epply  was  used  to  measure  solar  radiation.  A 
double  hemisphere,  high  precision  Epply  pyranometer  with  a  domed  Schott  RG8 
filter  as  its  outer  hemisphere  was  used  to  measure  near  infrared  radiation. 
The  reading  from  this  latter  pyranometer  was  corrected  according  to 
Drummond  and  Roche  (1965).  The  difference  between  solar  and  near  infrared 
radiation  was  assumed  equal  to  the  amount  of  visible  radiation.  Values 
of  visible  radiation  were  converted  from  energy  units  to  u  einsteins/cm^ 
sec  using  the  spectral  distribution  of  sun  and  sky  radiation  measured  by 
Federer  and  Tanner  (1966a). 

Leaf  Measurements 

As  indicated  in  the  theoretical  discussion,  the  leaf  area  and  leaf 
angle  distribution  as  functions  of  light  are  needed  in  the  radiation 
model.  These  variables  were  measured  using  an  improved  version  of  the 
grid  method  developed  by  Loomis  et  al.  (1969).  Corn  is  essentially  a 
two-dimensional  plant.  Thus,  it  is  possible  to  place  a  plant  in  a  natural 
position  against  a  large  grid.  In  this  case,  the  grid  was  made  up  of  1  dm 
squares.  Coordinates  were  noted  whenever  leaves  crossed  vertical  lines. 

The  variety  of  corn  (Cornell  M-3)  used  in  this  study  projected  its 
leaves  further  in  the  horizontal  than  the  vertical.  Thus  the  intersections 
with  the  vertical  lines  gave  a  more  detail  measure  of  the  leaf  shape  than 
intersections  with  the  horizontal  lines  as  used  by  Loomis  et^  al.  (1969). 
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Leaf  widths  were  measured  between  two  sets  of  coordinates.  Figure 
6  shows  a  typical  leaf  segment  against  part  of  the  grid.  The  leaf 
segment  between  (x-^.y-^)  and  was  approximated  by  a  straight  line. 

The  leaf  width  (W^)  was  measured  at  0.  The  area  (A)  of  this  segment  is 
expressed  in  d m2  as 

A  =  w1  •  /l  +  (y2  -  y^2  (33) 

The  angle  (IL)  of  the  leaf  segment  is  found  from 

IL  =  arctan  (y2  -  y^  (3U) 

In  figure  (6),  the  leaf  segment  crosses  a  horizontal  line.  Since 
a  knowledge  of  leaf  area  as  a  function  of  height  was  necessary,  A  must 
be  divided  into  A-^ (between  1-2  dm)  and  (between  2-3  dm).  These  are 
found  from 

A1  =  Wx  •  (x2  +  (2  -  y;L)2)1/2  (35) 

and  A2  =  W-l  •  ((l  -  x)2  +  (y2  -  2)2)l/2  (36) 

where  x  is  found  from  the  equation  of  a  straight  line, 

*2  “  xi 

x  =  X-.  +  -  (2  -  y  )  (37) 

1  y2  -  yx  1 

A  computer  program  was  written  to  do  the  above  calculations  for  each 
leaf  segment.  All  necessary  information  is  supplied  if  a  leaf  width 
is  provided  with  the  two  sets  of  coordinates  on  an  I.B.M.  card.  The  pro¬ 
gram  can  take  the  information  on  cards  and  calculate  the  leaf  area  and  leaf 
angle  distributions  as  functions  of  height  directly.  No  extensive  log 
for  each  plant  needed  to  be  made.  In  this  study  15-20  plants  were  measured 
for  each  plant  density  considered. 
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Results  and  Discussion 

The  purpose  of  the  visible  radiation  measurements  was  to  obtain  an 
average  irradiance  at  a  particular  level  in  the  crop  to  compare  with 
calculations  from  the  model.  An  average  value  (AVE)  is  defined  as: 

D 

/  R  d 

4  x  X 

AVE  =  - 9 -  (38) 

D 

/  dx 
0 

where  Rx  is  the  response  of  the  instrument  at  a  particular  distance 
along  the  sampling  track.  is  the  length  of  this  track.  Supposedly 
the  area  covered  by  the  width  of  the  sensing  head  times  is 
representative  of  the  crop  for  that  particular  level.  In  actual  fact 
the  above  integration  is  replaced  by  spot  readings  at  various  positions 
on  the  track.  Since  these  readings  are  at  regular  small  finite  intervals 
one  can  write ; 

D  N  N 

/  Rx  A  .E  Rx  *  A*  ,E  Rx 

AVE  =— -  ~  -Ai -  =  -ifi - 

A 

/  d  nax  N 

0 


where  N  is  the  total  number  of  observations . 

One  can  calculate  a  series  of  running  averages  for  increasing 
lengths  along  the  track  (Figure  7).  The  fact  that  the  averages  approach 
a  constant  value  as  x  -*■  D  assures  obtaining  a  representative  value  for 
the  crop. 

Measurements  of  the  plant  structure  are  needed  in  the  model.  Figures 
(8-12)  show  measurements  of  leaf  area  densities,  cumulative  leaf  area 
indexes  and  cumulative  leaf  angle  distributions  for  several  densities  of 
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Figure  8.  Cum.  L.A.I.  and  leaf  area  density  for  two  stands  of  corn  as 
functions  of  height  (1968  study). 

(A)  represents  cum.  L.A.I.  of  the  stand  with  total  L.A.I.  =  3. 
(o)  represents  cum.  L.A.I.  of  the  stand  with  total  L.A.I.  =  2. 
(•)  represents  leaf  area  density  (^S§.)  of  the  stand  with  total 

dm|  L.A.I.  =  3. 

(+)  represents  leaf  area  density  (jgfj  of  the  stand  with  total 
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LEAF  ANGLE  (IL) 


Figure  9.  Cumulative  leaf  angle  distribution  for  the  stand  of  L.A.I 
2.60.  The  leaf  area  is  divided  into  increments  of  1-6  dm 
7-9  dm,  10-13  dm  and  >  13  dm  (Study  made  in  1968). 
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Figure  10.  Cumulative  leaf  angle  distribution  for  the  stand  of  L.A.I. 

3.63.  The  leaf  area  is  divided  into  increments  of  1-7  dm, 
8-11  dm,  12-lH  dm  and  >  lU  dm  (Study  made  in  1968). 
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corn.  The  leaf  angle  distributions  are  of  particular  interest.  Leaf 
area  is  distributed  approximately  equally  over  the  angle  classes.  The 
distributions  tend  to  be  slightly  planophil  (more  horizontal  leaves  than 
erect)  especially  near  the  ground.  These  measurements  contrast  with  the 
spherical  distributions  Nichiporovitch  (i960)  found  for  corn. 

Figure  (13)  represents  calculations  using  both  the  de  Wit  and 
Duncan  models  for  a  leaf  area  index  of  3.63  and  a  sun  angle  of  52°.  The 
de  Wit  and  Duncan  models  give  approximately  the  same  answers  for  small 
values  of  S.  Reducing  S  from  0.202  to  O.lUO  made  little  difference  in  the 
de  Wit  model.  All  calculated  values  agree  reasonably  w^ll  with 
measurements. 

Figures  (lU-17)  represent  calculations  of  the  Duncan  Model  compared 
with  actual  measurements  in  several  stands  of  com  at  various  sun  angles. 
The  measurements  differ  from  one  another  not  only  for  sun  angle  but  also 
in  the  proportion  of  incoming  diffuse  or  sky  radiation.  While  attempts 
were  made  to  measure  only  when  the  sky  was  relatively  cloud  free,  especi¬ 
ally  near  the  sun,  perfectly  clear  skies  are  rather  rare  at  the  location 
of  the  field  site.  I  j 

1  ' 

The  1967,  data  were  taken  in  a  stand  of  corn  which  had  a  higher 
L.A.I.  and  in  general  was  more  uniform  than  the  1968  stands  of  com. 

Thus,  Figure  (lit)  shows  slightly  better  agreement  with  the 
model  than  the  1968  studies  (Figures  15,  l6). 

In  the  1968  measurements,  certain  trends  are  noticeable.  Measure- 

} 

ments  below  a  L.A.I.  of  2.  tend  to  be  higher  than  the  calculations. 

These  deviations  tend  to  increase  with  lower  sun  angles.  The  lowest 
angle  reported  shows  a  wide  deviation  from  the  model  but  these  low 
irracLance  levels  are  much  more  difficult  to  measure  exactly. 
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Figure  15.  Per  cent  transmission  as  a  function  of  CUM.  L.A.I.  at  various 
sun  angles  (Study  made  in  1968  with  L.A.I.  =  3.63  and  2.60). 
(•)  -  Duncan  model  for  corn  with  L.A.I.  =  3.63 

(O)  -  Measurements  for  corn  with  L.A.I.  =  3.63 

(o)  -  Duncan  model  for  corn  with  L.A.I.  =  2.60 

(A)  -  Measurements  for  corn  with  L.A.I.  =  2.60 
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Figure  l6.  Per  cent  transmission  as  a  function  of  CUM.  L.A.I.  at  various 
sun  angles  (Study  made  in  1968  with  L.A.I.  =  3.13  and  2.38). 
(•)  Duncan  model  for  corn  with  L.A.I.  =  3.13 

(o)  Measurements  for  corn  with  L.A.I.  =  3.13 

(+)  Duncan  model  for  corn  with  L.A.I.  =  2.38 

(A)  Measurements  for  corn  with  L.A.I.  =  2.38 
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Die  Duncan  model  assumes  random  distribution  of  leaf  area.  Com 
planted  in  a  15"  hexagonal  pattern  with  plants  oriented  randomly  about 
their  stems  should  give  a  good  approximation  to  this  random  distri¬ 
bution.  Any  point  in  the  canopy  could  have  leaves  from  four  to  five 
plants  in  its  vicinity.  Another  assumption  is  independence  between 
leaf  increments.  If  a  plant  is  missing  from  the  above  pattern  which 
invariably  happens  in  large  plots  of  machine  planted  com,  a  certain 
amount  of  dependence  between  levels  is  introduced.  Diis  would  allow 
more  light  into  the  lower  levels  than  the  model  would  predict  and  would 
give  a  larger  error  for  com  with  a  lower  L.A.I. 

Figure  (17)  is  particularly  interesting  as  it  represents  measure¬ 
ments  from  a  completely  overcast  sky.  Die  agreement  of  measurements 
with  the  model  is  particularly  encouraging  for  this  rather  extreme  type 
of  cloud  cover.  However,  the  effects  of  reflection  from  cloud  banks  for 
partially  overcast  skies  have  not  been  studied  here  and  could  lead  to 
serious  errors  for  partially  cloudy  days  . 

Die  scattering  theory  is  tested  in  figure  (l8 ) .  Where  the  propor¬ 
tion  of  reflected  radiation  is  plotted  as  a  function  of  sun  angle  for 
both  calculated  and  measured  values.  Diere  is  reasonably  good  agreement 
between  the  model  and  calculations  at  sun  angles  from  55°  to  65°.  Diere 
is  a  tendency  for  the  model  to  overestimate  at  the  highest  sun  angles 
(6O-650).  As  sun  angle  decrease  below  50°,  the  model  underestimates  the 
measured  values  in  increasing  amounts.  This  deviation  is  probably  due 
to  specular  ’eflection  from  the  leaf  surfaces.  In  the  model  we  assume 
all  light  is  reflected  diffusely.  Because  the  amounts  of  radiation 
involved  are  quite  small,  these  discrepencies  shouldn't  influence  a 
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photosynthesis  model.  However  the  scattering  theory  introduced  in  this 
study  is  not  good  enough  when  the  measure  of  reflected  radiation  is  the 
prime  concern  such  as  in  remote  sensing  situations. 


Conclusions 


The  radiation  model  based  on  a  random  distribution  of  leaf  area 
successfully  predicts  the  penetration  of  visible  radiation  into  com. 

There  is  a  tendency  to  overestimate  amounts  of  radiation  at  the  bottom 
of  the  stands  at  the  lower  leaf  area  indexes.  However,  this  was  not 
considered  a  serious  error.  In  fact,  one  could  argue  that  if  the  leaf 
area  surrounding  gaps  in  the  planting  pattern  is  small  compared  with 
the  total  leaf  area,  the  model  will  describe  the  light  regime  better  than 
the  measurements  sis  fsir  as  calculations  for  photosynthesis  sire  concerned. 

The  bsisic  equations  relating  percent  transmission  to  cum.  leaf 
area  index  from  both  the  de  Wit  and  Duncan  models  were  tested  in  this 
study.  They  gave  essentially  the  same  answers  when  S,  the  amount  of 
leaf  sirea  per  increment,  less  than  0.2  was  used.  The  Duncan  equation  was 
favored  because  it  wsis  independent  of  S  and  thus  less  empirical.  Also 
final  answers  could  be  obtained  with  fewer  calculations  as  the  leaf  area 
could  be  divided  into  larger  increments. 

In  general,  the  agreement  between  equations  and  measurements  was 
sufficiently  good  that  a  general  photosynthesis  model  was  constructed 
based  on  the  development  outlined  above.  Variolas  parts  of  this  model  will 
be  considered  in  subsequent  chapters. 


Chapter  IT 


A  LEAF  PHOTOSYNTHESIS  MODEL 

In  many  instances ,  one  can  best  treat  a  crop  in  terms  of  an  array 
of  leaves.  Thus,  it  is  particularly  useful  to  first  mathematically 
model  the  responses  of  individual  leaves.  This  process  incorporates  a 
large  amount  of  information  into  the  concise  form  of  mathematical  equations. 
The  information  is  then  readily  accessible  when  one  integrates  over  the 
total  leaf  area  in  a  crop  response  model. 

In  this  study  a  mathematical  treatment  is  developed  to  describe 
net  photosynthesis  of  leaves  as  a  function  of  temperature,  CO2  concen¬ 
tration  at  the  leaf  surface  and  light  flux  density.  At  this  stage,  the 
model  is  highly  specific  to  com  and  other  tropical  grasses.  The  model 
is  tested  against  experimental  leaf  chamber  measurements.  The  measure¬ 
ments  were  made  on  attached  leaf  blades  of  field  grown  com. 


Theoretical  Considerations 


Representing  gaseous  exchange  characteristics  of  a  leaf  by  a  network 
of  resistances  is  becoming  increasingly  common.  (For  example.  Moss 
(1966),  Lake  (1967)  and  Waggoner  (1969)).  The  representation  by  Lake 
is  the  most  detailed,  and  with  some  modification,  is  presented  in  Figure 
(19).  In  this  figure,  rg  is  the  stomatal  resistance,  r^  is  the  resis¬ 
tance  of  the  cell  walls,  rm,  rwl,  and  r^  are  internal  cell  resistances 
and  rc  is  a  carboxylation  resistance  to  be  defined  below. 

A  com  leaf  under  a  visible  radiation  load  releases  very  little  or 
no  CO^  into  a  C02  free  environment  even  when  the  stomates  are  open  (Moss 
(1966)).  It  is  generally  agreed  that  r_ ,  r^  and  r  are  of  the  same  order 
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of  magnitude,  (Waggoner  (1969)).  Therefore  for  com,  rw2  is  much  greater 
than  rwl  and  Figure  (19)  can  be  simplified  to  Figure  (20 ). 

Following  Chartier  ( 1966b ),( 1969 )  and  using  Figure  (20),  net  photo¬ 
synthesis  (N)  and  gross  photosynthesis  (P)  can  be  expressed  as 


c  -  cx 

r  +  r 

(39) 

s  m 

C  -  C, 

1  +  R 

(to) 

rs  +  rm 

where  R  is  respiration  and  C-^  is  the  concentration  of  CO2  at  the  chloro- 
plasts.  Rabinowitch  (1951)  derived  an  expression  relating  Cj,  light 
irradiance (i)  and  P  by  assuming  the  following  equations; 


P  =  k-jA]  Cx 

(1*1) 

P  =  k2[AC02l  [H] 

(1*2) 

Site  of 
photosynthesis 


Site  of 
respiration 

Figure  19.  Hie  resistance  network  for  a  leaf:  a  modification  of  the  network 
by  Lake  (1967). 
rs  -  stomatal  resistance 
ri  -  cell  wall  resistance 
rm,  iv!,  rW2  "  internal  cell  resistances 
rc  -  carboxylation  resistances. 


Figure  20.  A  simplified  resistance  network  for  a  leaf 
N  -  net  photosynthesis 
P  -  gross  photosynthesis 
R  -  respiration 


•  '  •  •;  . 


1*9 


and 


[H]  =  k3I  . 


(1*3) 


[A],  [ACO2]  and  [H]  are  concentrations  of  ribulose  diphosphate,  phos- 
phoglyceric  acid  and  triphosphopyridine  nucleotide  respectively,  k^,  kg 
and  k^  ar^:  proportionality  constants.  Combining  equations  (1*1-1*3)  one 
obtains ; 

ol 


P  = 


1  +  I 


where  a  =  k„  k.,  A0 


and  r  = 
c 


2  3 
1 

Aokl 


(1*1*) 

(1*5) 

(1*6) 


When  is  eliminated  from  equation  (ho)  using  (1*1*),  the  following 
quadratic  is.  obtained: 

(rs  +  r^P2  -  P[al(rs  +  rm  +  rc)  +  C  +  R(rs  +  rm)]  +  Cal  +  Ral(rg  +  rm)  =  0 

(1*7) 

The  stomatal  resistance  under  conditions  of  no  water  stress  often  appears 
as  a  hyperbolic  function  of  I  of  the  form; 


rs  =  Y  + 


0 


I  +  I' 


(U8) 


where  y,  0  and  I’  are  constants  (Kuiper  (l96l)).  Its  incorporation  here 
is  an  improvement  on  Chartier's  (1966b)  treatment  where  he  assumed  rg  was 
a  constant.  In  equation  ( 1+8)  the  differential  of  P  with  respect  to  I 
at  P  and  I  equal  to  zero  is  found  as 


3P 

31 


=  a 


(1*9) 


P  =  0 
1  =  0 


1  1 1  I  i  ■  ■  ! 1 "  I  ■  J'l  >f  l>U  *■■■ 
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When  I  is  large,  P  -*■  P..  and  r  ■*■  y 
Equation  (4?)  reduces  to 


C  +  R(y  +  r 

P  =  _ 2_ 

fm 

Y  +  rm  +  rc 

Let  NM  =  PM  -  R 

Then  H  .  C  1  B>Y  1  r->  I  B(V  1  1  r°> 

M 

Y  +  r  +  r 
m  c 


R  r 


or 


N  = 
M 


Y  +  r  +  r  Y  +  r„  +  r„ 
me  me 


(.50) 

(51) 

(52) 

(53) 


Hie  second  term  on  the  right  in  (53)  is  small  and  a  good  approximation  of 
Nm  is  given  by; 


nm  = 


y  +  r  +  r 
1  m  c 


(54) 


Finally  by  adding  rQ  to  both  sides  of  equation  (47)  one  obtains 


rc  =  {rs  +  rm  +  rc)  ~ 


aI(rs  +  rm  +  rc>  +  c  +  R(rs  +  rm} 


Cl a  t  Ral(r  +  rj 
s  m 


(55) 


Thus  rc,  and  Y  +  rm  +  rc  can  be  found  experimentally  from  a  series  of 
measurements  of  net  photosynthesis  at  various  light  flux  densities  at  a 
known  concentration  of  CO^.  The  procedure  will  be  considered  in  more 
detail  below. 

In  this  treatment,  photosynthesis  (P)  is  assumed  independent  of 


temperature.  However,  temperature  effects  on  N  are  introduced  via  the 


•  •*  «*«-*#*  n*  „ 
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respiration  term  (r).  Following  Waggoner  (1969),  the  respiration  rate 
can  be  expressed  as 

R  =  Rx  EXP [9000  ln(Q)  (—  -  |)  ]  (56) 

where  Rx  is  the  respiration  rate  at  303°K.  Ihe  temperature  T  is  in  °K. 

Q  is  the  rate  of  increase  of  R  for  each  10°  rise  in  temperature  and  is 
assumed  equal  to  2.  Note  that  "photorespiration"  has  been  ignored 
here  which  again  limits  this  development  to  com  and  other  efficient 
tropical  grasses  (Moss  (l 966)). 

In  summary,  net  photosynthesis  (N)  can  now  be  expressed  as  the 
following  function  of  light  irradiance  (i),  CO concentration  at  the  leaf 
surface  (c)  and  temperature  (T); 


..  -Bx  -  Ax  -  '•V* 

N  =  - —  -  R 

2A 

x 


where  R  =  RY  EXP[9000  InQ  (—  -  -) 

303  T 


/ 


rs  +  rm 


-(al(rs  +  rm  +  rc)  +  C  +  R(rs  +  rm) ) 


(57) 

(58) 

(59) 

(60) 


and  Cx  =  Cal  +  Ral(rs  +  rm) 


(61) 


Materials  and  Methods 


Leaf  chamber  measurements  were  made  by  Dr.  R.  B.  Musgrave  of  the 
Cornell  University  Agronomy  Department.  The  leaf  chamber  has  been 
described  in  detail  by  Heikel  (1968).  In  brief,  the  chamber  is  of  the 
semi-open  type.  That  is,  air  enters  an  inlet  port  to  the  leaf  chamber 
and  out  via  an  outlet  port  as  in  an  open  system.  While  in  the  leaf 
chamber,  air  is  rapidly  recirculated  through  the  chamber  and  a  set  of 
cooling  coils.  This  recirculation  keeps  the  air  at  a  high  turbulence 
level  and  provides  control  of  the  leaf  temperature.  Also  the  air  inside 
the  chamber,  because  of  the  rapid  mixing,  is  at  a  uniform  concentration 
of  CC>2.  Net  photosynthesis  is  calculated  from  me£isurements  of  the  flow 
rate  at  the  inlet  and/or  outlet  port  and  the  CO^  concentration  difference 
between  the  inlet  and  outlet  ports.  The  C02  concentration  inside  the 
chamber  can  be  varied  by  adjusting  the  flowrate  into  the  chamber.  Input 
air  is  drawn  from  a  large  mixing  tank  to  avoid  the  usual  fluctuations  of 
COg  concentration  in  the  ambient  air  under  field  conditions.  The  leaf 
chamber  is  large  enough  to  enclose  one  whole  leaf-blade.  It  is  also 
quite  mobile  so  that  attached  leaves  of  field  corn  can  easily  be  measured. 
Leaf  temperature  is  controlled  to  a  maximum  variation  of  ±2°C(Heikel  and 
Musgrave  (1969)).  In  this  study,  the  radiation  load  was  supplied  by  a 
bank  of  incadescent  lamps  (Sylvania,  #300PAR56/2NSP) .  These  particular 
lamps  transmit  60  percent  of  the  infrared  radiation  through  the  back  of 
the  lamp.  This  considerably  reduces  the  heat  load  on  the  leaf.  Light 
irradiance  levels  on  the  leaf  were  changed  by  placing  layers  of  cheese¬ 
cloth  between  the  lamp  and  the  leaf  chamber.  The  irradiance  in  the 
chamber  under  various  numbers  of  cheesecloth  layers  was  measured  by  an 


uncorrected  Weston  sun-light  meter.  The  readings  were  converted  to 
U  einsteins/cm^sec  using  conversion  factors  of  Gaastra  (1959)* 

Measurements  were  made  at  the  Microclimate  Investigations  Experi¬ 
mental  Site  in  Ellis  Hollow  near  Ithaca,  New  York  during  July  and  August, 
1968.  The  variety  measured  was  Cornell  M-3.  The  measurements  were  made 
in  the  same  field  as,  and  in  cooperation  with,  studies  of  natural  exchanges 
of  CO^,  heat  and  water  vapor  between  a  com  crop  and  the  environment 
conducted  by  Dr.  E.  R.  Lemon  and  associates,  individual  leaf  measure¬ 
ments  of  net  photosynthesis  at  full  sunlight  intensity  were  made  on 
various  leaves  of  several  plants.  More  important  to  this  study  were 
measurements  of  net  photosynthesis  as  a  function  of  light  intensity  at  a 
constant  concentration  of  CO2  on  two  separate  leaves.  Net  photosynthesis 
as  a  function  of  C02  concentration  at  a  constant  light  flux  density 
(0.230  U einsteins/cm2sec)  was  measured  on  one  leaf.  Respiration  rates 
at  30°C  were  also  measured  in  the  leaf  chamber. 

Field  measurements  of  stomatal  resistance  as  a  function  of  light 
flux  density  were  rade  independently  of  the  leaf  chamber  measurements  by 
R.  W.  Shawcroft"*".  A  Van  Bavel  porometer  (Van  Bavel  et_  al.  (1965))  was 
used  to  measure  stomatal  resistance  to  water  vapor  diffusion  of  leaves 
in  their  natural  positions  in  the  field.  The  resistance  to  diffusion  of 
CO^  was  obtained  by  multiplying  the  resistances  to  water  vapor  by  the 
ratio  of  the  diffusion  coefficients  for  water  vapor  and  for  CC>2  in  air. 
Measurements  were  made  at  various  times  on  a  clear  day  on  both  shaded 
and  unshaded  leaves  to  obtain  stomatal  resistances  at  8  large  range  of 
light  flux  densities.  Resistances  of  the  upper  and  lower  surface  of  the 
leaf  were  added  in  parallel  to  give  the  total  resistance  of  the  lei_f. 

^Unpublished  data  to  be  submitted  in  partial  fulfillment  of  the 
Ph.D.  requirements  at  Cornell  University. 


Immediately  after  making  a  stomata!  resistance  measurement  the  light 
sensor  described  in  Chapter  I  was  placed  parallel  and  close  to  the  leaf 
plane  to  obtain  the  light  flux  density  at  the  leaf  surface.  A  similar 
procedure  was  used  by  Turner  (1969).  The  data  were  fitted  by  regression 
to  the  hyperbola  described  by  equation  (1+8)  to  obtain  values  of  y ,  3  and 
I'.  The  measurements  were  made  when  soil  moisture  was  not  limiting. 

This  theory  is  extended  in  Chapter  IV  to  conditions  of  water  stress  by 
increasing  the  value  of  y. 
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Results  and  Discus s i ons 


In  equations  (57-61),  unknown  leaf  parameters  are  it,,  rgj  r  r 
and  «.  Bx  is  the  leaf  respiration  at  303°K.  From  measuresmnts  in  a  dark 
leaf  ohamber,  Rx  was  found  to  be  5.0  mg/dm2  hr.  Die  stomatal  resistance 
(rs)  was  related  to  the  light  flux  density  by  equation  (16) 


rs  =  7  + 


Itl' 


From  the  field  measurements  by  R.  w.  Shaweroft,  Y,  6  and  I'  ware  found 
to  be  1.146  sec/cm,  0.036  seAm/u  einstein  and  0.0015  u  einsteinsW  sec 
respectively  for  the  com  variety,  Cornell  M-3,  used  in  this  study. 

Fhotosynthesis  as  a  function  of  light  irradiance  is  used  to  determine 
r"  ™'1  rc'  FleUre  l21>  sh°"s  measurements  of  H  at  various  light  irrad- 
iances.  A  curve  is  drawn  by  hand  through  the  points.  From  this  curve 
estimates  of  a  and  are  determined  graphically. 

Let  r  =  r  +  t  +  -r 

s  m  rc  (62) 

From  (L8),  (5I+)  and  (62) 

_  _c  _ e_ 

»M  i*r  (63> 

Knowing  »M  from  Figure  (21),  r  is  determined  as  a  function  of  I.  Com- 
bining  (62)  and  (59)>  one  obtains 


rc  "  r" 


“Ir  +  C  +  R(rs  +  r  )  Clot  +  Ral(r  +  rJ 


Hius  rc  can  be  evaluated  for  a  given  rm  using  values  from  the  curve  of 


1.0 


0  y  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.9  0.9 

I  (fraction  of  0.2304  jj  einsteins/cm2.sec) 


Figure  21.  Net  photosynthesis  (N)  as  a  function  of  light  flux  density  (i) 
at  a  CC>2  concentration  of  260  ppm. 
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Figure  (21).  Finally  rffi  is  evaluated  up  g 

rm  =  r  -  rs  -  rc  (65) 

Equations  (64)  and  (65)  are  solved  simultaneously  for  r  and  r  .  In 

m  c 

practice  rc  and  rm  are  not  constants  with  I.  This  represents  a  discrep- 

ency  between  theory  and  measurements  since  both  rc  and  rm  are  defined 

as  constants.  Ten  values  of  r  and  r  are  evaluated  at  light  irradiances 

c  n 

of  0.1,  0.2,  0.3  .  •  .  of  0.230  y  einsteins/cm^sec.  (0.230  y  ej.nsteins/ 
cm^sec  corresponds  to  10,000  ft.  candles).  These  values  are  averaged 
and  used  in  equations  (57-6l)  to  generate  theoretical  light  and  C02 

p 

curves.  By  using  a  value  of  of  47.4  mg/dm  hr  which  is  a  reasonable 
estimate  from  Figure  (2l),  and  average  values  of  rm  and  rc  the  theory 
tends  to  overestimate  N  at  intermediate  light  flux  densities  (a  maximum 
of  16$)  and  underestimate  N  at  high  flux  densities  (a  maximum  of  h%) .  These 
differences  were  reduced  to  Q%  and  1 %  respectively  by  adjusting  NM  from 
47.4  to  49.4  mg/dm  hr.  Corresponding  values  of  ?c  and  rm  for  NM  equal 
to  49.4  mg/dm2  hr  are  0.20  and  I.65  sec/cm  respectively.  This  adjustment 
also  improved  agreement  between  the  theoretical  and  measured  C02  curve. 

These  theoretical  curves  with  measurements  are  shown  in  Figures  (22  )  and 
(23). 

It  should  be  pointed  out  that  measurements  in  Figure  (23)  were  made 

O 

at  a  light  irradiance  of  0.230  y  einsteins/cm  sec.  The  leaf  used  in 
Figure  (23)  had  a  higher  net  photosynthesis  rate  (N  =  56.5  mg/dn2  hr)  at 

p 

260  ppm  C02  than  the  two  leaves  used  in  Figure  (22  )  (N  =  46.5  mg/dm  hr). 

The  theoretical  values  in  Figure  (23)  were  obtained  by  multiplying 
equation  (57)  by  where  Pr  =  56.5/47.0.  When  this  leaf  model  is  used 


Figure  22.  Net  photosynthesis  (N)  as  a  function  of  light  flux  density  (i) 
at  a  CO2  concentration  of  260  ppm. 
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Figure  23.  Wet  photosynthesis  (N)  as  a  function  of  C02  concentration  at 
a  light  flux  density  of  0.230  p  einsteins/cnr/sec. 


in  the  simulation  of  the  crop  response  is  defined  as 


P 


r 


^260 

1+6.5 


N26o  is  the  value  of  net  photosynthesis  measured  at  260  ppm  C0o  at 
0.230  u  einsteins/cm^sec  averaged  over  a  representative  number  of  leaves 


in  the  crop.  Note  that  if  is  measured  as  a  function  of  depth  into 

the  crop,  this  function  can  easily  be  incorporated  into  the  crop  model. 


The  measurements  made  in  this  study  showed  no  significant  variation  of 
N  with  depth.  P2g0  was  found  to  be  1+7.2  mg/dm2hr  for  the  field  corn 


under  consideration  in  this  study. 


-  »waa 
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Conclusions 

A  mathematical  model  of  net  photosynthesis  of  an  individual  corn 
leaf  has  been  developed  in  this  study.  The  model  is  tested  using 
experimental  leaf  chamber  measurements  of  net  photosynthesis  as  functions 
of  C0o  concentrations  and  light  flux  densities.  By  forcing  r  ,  the 
car-boxy  lation  resistance,  and  rm>  the  mesophyll  resistance,  to  be  con¬ 
stants,  the  theory  tended  to  slightly  overestimate 'photosynthesis  rates 
at  intermediate  light  flux  densities.  However  the  agreement  between 
theory  and  measurements  was  considered  satisfactory  (a  maximum  discrepency 
of  8$)  considering  the  rather  simple  assumptions  used  in  the  theoretical 


model . 


Chapter  III 


THE  BOUNDARY  LAYER  RESISTANCE  OF  PLANT  LEAVES 

The  concept  of  a  boundary  layer  was  first  conceived  by  Prandtl  who 
divided  flow  around  a  body  into  two  parts ;  a  thin  layer  near  the  surface 
of  the  body  where  friction  or  shear  plays  an  important  role  and  the  remain¬ 
ing  part  of  the  flow  where  friction  may  be  neglected.  Boundary 
layer  theory  which  was  developed  around  this  concept  has  proved  extremely 
useful  in  engineering  applications.  It  is  not  surprising  therefore  that 
these  ideas  have  been  applied  to  flow  of  air  over  leaves.  A  leaf  in  its 
natural  environment  exchanges  heat,  COg,  Og  and  water  vapor  with  what 
is  normally  a  highly  turbulent  air  stream.  The  resistance  to  diffusion 
encountered  in  the  leaf  boundary  layer  is  relatively  large  compared  to 
eddy  diffusion  in  the  turbulent  air  stream.  For  this  reason  it  is  an 
important  consideration  in  the  study  of  the  energy  balance  of  single  leaves, 
net  photosynthesis  and  the  water  use  of  plants. 

The  concept  of  diffusion  of  mass  and  heat  across  boundary  layers 
is  expressed  in  terms  of  a  resistance  to  diffusion,  r&.  Resistances  to 
the  diffusion  of  heat  (ra^),  water  vapor  (ray),  and  CO2  (rac)  are  defined  as 


rah  =  PCp 


pCp 

raw  =  - 

Y' 


LE 


and  rac 


(66) 

(67) 

(68) 


62 


63 


H,  LE  and  N  are  respectively  fluxes  of  sensible  heat,  latent  heat  and  net 
photosynthesis  per  unit  area,  e^,  and  are  vapor  pressure, 
temperature  and  CCL  concentration  at  the  leaf  surface:  e„,  TD  and  CL  are 
the  corresponding  values  of  the  air  stream;  p,  Cp  and  y'  are  the  density 
of  air,  the  specific  heat  of  air  and  psychrometric  constant  respectively. 
This  concept  of  resistance  lends  itself  to  the  understanding  of  exchange 
processes  from  a  leaf.  For  example ,  rac  can  be  added  to  the  left  of  rs 
in  the  series  of  resistances  of  Figure  ( 20)  in  Chapter  II.  The  leaf 
itself  has  two  surfaces  ,  and  ra  for  a  given  surface  should  be  added 
in  series  with  the  stomatal  resistance  for  that  surface.  However,  measure¬ 
ments  by  Shawcroft  described  in  Chapter  II  indicated  that  the  stomatal 
resistance  for  the  upper  and  lower  leaf  surfaces  were  approximately  equal. 
Thus ,  only  a  small  error  is  introduced  by  adding  the  total  leaf  stomatal 
resistance  and  the  total  leaf  boundary  layer  resistance  in  series.  In  a 
leaf  chamber  rac  is  considered  zero  because  of  the  highly  turbulent 
flow  in  the  chamber.  In  the  natural  environment  this  is  not  necessarily 
true . 

In  order  to  extrapolate  leaf  chamber  measurements  to  the  natural 
field  environment,  the  size  of  the  leaf  boundary  layer  resistance  under 
natural  conditions  must  be  approximated.  From  dimensional  analysis,  the 
boundary  layer  resistance  is  a  function  of  a  length  parameter  (in  this 
case  the  width  of  the  leaf),  the  velocity  of  the  free  air  stream  in  the 
vicinity  of  the  leaf  and  in  the  case  where  free  convection  is  important, 
to  the  temperature  difference  between  the  leaf  and  the  air  stream.  To 
obtain  these  functions  experimentally  one  is  confronted  with  two  alterna¬ 
tives.  The  experimenter  can  work  in  a  wind  tunnel,  usually  with  artificial 
leaves,  and  determine  precise  relationships  (For  example  Thom  1968, 


6k 


Parkhurst  et_  al_.  (1968)).  However,  since  the  level  of  free  stream  tur¬ 
bulence  is  very  low  in  wind  tunnels,  the  experimenter  has  difficulty 
relating  his  measurements  to  field  conditions.  On  the  other  hand, 
working  in  the  field  on  intact  leaves  also  presents  problems.  Here  the 
measurements  are  necessarily  less  precise  since  working  conditions  are 
far  from  ideal.  Also  unsteady  state  conditions  in  the  field  make  it 
difficult  to  relate  measurements  to  the  wind  velocity  and  leaf  to  air 
temperature  gradients. 

Two  studies  which  were  made  on  attached  leaves  in  the  field  are 
Hunt  et  al.  (1968),  and  Kanemasu  et_  al.  (1969).  An  energy  balance  method 
was  used  in  both  studies.  (Hunt  et  al.  (1968),  also  used  a  transient 
method  which  will  be  discussed  in  detail  below).  Under  steady  state 
conditions  the  energy  balance  equation  of  a  leaf  is 


PCp  E(Tl)  -  Ea 
Y'  raw  +  rs 


(69) 


where  is  the  net  radiation  load  on  the  leaf,  rg  is  the  stomatal  re¬ 
sistance  to  water  vapor  flow  and  E(T^)  is  the  saturated  vapor  pressure 
as  a  function  of  leaf  temperature.  If  one  assumes  ra„.  =  ra^  and  measures 
all  other  components  [Rj^  is  measured  with  a  net  radiometer  positioned 
above  and  below  the  leaf;  rs  is  measured  with  a  leaf  porometer  (eg.  Van 
Bavel  et  al . ,  1965)],  equation  (69)  can  be  solved  for  ra.  However,  a 
large  number  of  components  have  to  be  measured  before  ra  can  be  cal¬ 
culated.  This  necessarily  reduces  the  accuracy  of  the  method.  Kane¬ 
masu  et  al.  (1969)  reduced  the  number  of  variables  by  measuring  the  latent 
heat  term  directly.  All  but  the  leaf  under  study  were  stripped  from  a 
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bean  plant  growing  in  a  pot.  The  pot  was  weighed  at  intervals  to  deter¬ 
mine  the  latent  heat. 

Most  of  the  wind  tunnel  work  to  date  has  been  with  artificial  leaves. 

For  example,  Parkhurst  et  al.  (1968)  made  artificial  leaves- out  of  two 

layers  of  thin  copper  sheeting  with  heating  wire  between  the  copper 

sheeting.  By  passing  a  constant  known  current  through  the  heating  wire 

and  measuring  T^  and  T&  with  thermocouples,  equation  (66)  was  used  to 

determine  ra^.  Thom  (1968)  used  a  similar  technique  to  measure  ra^. 

In  addition  he  glued  blotting  paper  to  his  metal  leaves  and  soaked  them 

with  various  liquids.  By  attaching  the  metal  leaf  to  an  arm  of  a  weighing 

balance,  evaporation  from  the  leaf  was  estimated.  Thus  Thom  was  able 

to  measure  ra  as  well  as  ra.  . 

w  n 

There  are  some  properties  of  leaves  such  as  leaf  hairs,  curvature 
and  flexibility  which  are  difficult  to  simulate  with  artificial  models. 
Pearman  (1965)  used  a  transient  method  to  measure  rah  on  real  leaves 
under  controlled  conditions .  A  leaf  was  heated  above  the  ambient  air 
temperature  with  microwaves.  Transient  times  were  recorded  from  the 
time  the  microwave  source  was  shut  off  until  leaf  temperature  equilibrium 
was  re-established.  Natural  and  forced  convection  was  studied. 

In  general  then,  there  are  few  data  on  boundary  layer  resistances 
of  real  leaves.  Because  of  this  and  a  need  for  specific  data  on  com 
(Zea  mays)  leaves,  a  transient  method  was  developed  to  measure  ra^  on 
excised  leaves  placed  in  a  wind  tunnel.  To  understand  this  transient 
method  some  theoretical  considerations  are  needed. 
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Theoretical  Considerations 


Consider  a  leaf  under  a  radiation  load  with  a  high  shortwave 
component.  If  the  radiation  load  changes,  the  leaf  temperature  will 
adjust  to  a  new  level.  The  final  equilibrium  state  can  be  described  by; 


h  pCp  .  .  pCP  e(Tf)  -  e 

R_  -  2ec  T^4  -  — -  (Tf  -  Ta) - ....  1  a 


nc 


rav 


Y'  ra^  +  rs 

where  ^  is  the  total  incoming  radiation,  Tj,  is  the  final  leaf  tempera¬ 
ture,  e  is  the  thermal  emissivity  of  the  leaf  (taken  to  be  1.0  in  this 
study)  and  o  is  the  Boltzman  constant.  The  unsteady  state  can  be 
described  by 

pCp  (e(TL)  -  ea) 


=  0.0 


(70) 


i  pCP 

Rn  -  2ca  Tl>  -  (Tl  -  Ta)  - 

s  rah  Y  raw  +  rs 


=  M  s  dTL  (7i) 
A  dt 


T^  is  the  leaf  temperature,  M/A  is  the  mass  per  unit  area  of  the  leaf,  s 
is  the  specific  heat  of  the  leaf  (taken  to  be  0.9) and  t  is  time.  Com¬ 
bining  equations  (70)  and  (7l)  one  obtains: 


P  Cp 

St  {tl  -  Tf }  + 


PCp 

y' 


e(TL)  -  e(Tf) 
(ra„-  +  rs) 


+  2ea(Tj} 


M  dTL 
A  S  dt 


(72) 


at  this  point,  two  simplifying  assumptions  are  made. 

First;  e(T)  =  bT  +  d  (73) 

where  b  and  d  are  constants.  This  is  a  reasonable  assumption  between 
20°  and  35°C.  In  this  range,  b  is  equal  to  2.08  mb/°C;  d  is  not  used  in 
the  calculations.  Also  one  assumes; 


T^  -  Tfk  =  >+T3  (Tl  -  Tf) 


(71*) 


where  T  is  the  average  of  T0  and  T^.  where  T0  is  the  initial  temperature. 


i  iimir  axiitUfifaM 


'Hie  percentage  error  of  this  assumption  can  be  expressed  as; 


Percentage  Error  = 


/  UeaT3(TL  -  Tf)  dt 


/  ea(T  **  -  TfS  dt 
0  L 


x  100. 


Equation  (75)  can  be  evaluated  by  assuming  an  exponential  increase  of 
decrease  of  leaf  temperature  with  time.  The  percentage  error  has  a  maxi¬ 
mum  value  of  k%  for  the  temperature  differences  considered  in  this  study. 
Using  the  above  simplifications ,  equation  (72)  becomes; 


(Tl  -  Tf)  (2eoT3  +  - - - +  —  )  =  -  !s 

y'  (raw  +  rs)  ra^  A 


d(TL  -  Tf) 


By  integrating  equation  (76)  and  solving  for  rah,  one  obtains; 


ra*  =  pCp/[—  s  i£(il  -  (8eoT3  ♦  - 5 - )] 

A  t1/2  y'  (ra^  +  rs) 


where  t^  is  the  time  required  for  the  leaf  to  reach  the  midway  point 
between  the  initial  and  final  leaf  temperatures.  By  assuming  ra„-  is 
equal  to  rah,  equation  (77)  can  easily  be  solved  by  iteration  for  ra^. 
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Materials  and  Methods 

A  com  leaf  was  cut  from  a  plant  under  water,  placed  in  a  beaker 
so  that  the  cut  end  was  under  water  and  inserted  into  the  working  section 
of  the  wind  tunnel  through  a  hole  cut  in  the  floor.  The  rim  of  the  beaker 
was  flush  with  the  floor  of  the  wind  tunnel.  A  leaf  was  arranged  to 
protrude  at  an  angle  from  the  floor  in  a  more  or  less  natural  manner  or 
was  placed  in  a  leaf  holder  so  that  a  large  portion  of  the  leaf  was  held 
at  zero  angle  of  attack  to  the  wind  velocity  across  the  width  of  the 
wind  tunnel.  A  third  treatment  consisted  of  placing  two  leaves  upwind 
of  the  leaf  in  the  natural  position.  This  was  done  to  study  the 
effect  of  free  stream  turbulence  generated  by  the  two  upwind  leaves  on  the 
boundary  layer  resistance  of  the  third  leaf.  Turbulent  intensities  were 
measured  by  a  heated  thermocouple  annemometer.  (Hastings,  Model  Am-62X). 

A  bank  of  incadescent  lamps  immersed  in  a  plexiglass  water  bath  was 
placed  on  the  ceiling  of  the  working  section.  The  ceiling  and  sides  of 
the  working  section  were  made  from  clear  plexiglass  as  well.  Thus  a 
shortwave  radiation  load  equivalent  to  full  sunlight  intensity  could  be 
placed  on  the  leaf.  One  could  change  the  radiation  intensity  by  changing 
the  setting  of  a  large  rheostat  which  was  connected  in  series  with  the 
lights . 

Very  fine  copper-constantan  thermocouples  (.005  inches  in  diameter) 
were  used  to  measure  the  leaf  to  air  temperature  difference.  Larger 
copper-constantan  thermocouples  referenced  to  an  ice  bath  were  used  to 
measure  the  air  temperature  at  each  end  of  the  working  section  of  the 
windtunnel.  All  temperatures  were  continuously  recorded  on  strip  chart 
recorders.  By  changing  the  radiation  load  on  the  leaf  and  recording  the 
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leaf  to  air  temperature  difference  continuously,  could  be 
evaluated.  The  air  temperature  remained  reasonably  constant  for  small 
changes  in  the  radiation  loads.  Ibis  procedure  was  repeated  six  to  ten 
times  at  various  radiation  loads  for  each  wind  speed.  Each  treatment 
was  repeated  for  three  to  four  leaves. 

Equation  (77)  has  a  stomatal  resistance  tea,  rs,  which  his  not 
been  discussed  as  yet.  leaves  were  placed  in  the  wind  tunnel  prior  to 
sunset  and  left  in  the  dark  for  several  hours.  This  was  done  to  close 
the  stomates  prior  to  the  start  of  the  experiment.  The  experiment  was 
conducted  at  night  since  stomates  of  com  and  many  other  agricultural 
crops  once  closed  do  not  open  readily  at  night  under  a  radiation  load. 
Preliminary  measurements  made  with  a  porometer  indicated  that  ^  values 
were  of  the  order  of  lU.O  sec/cm  for  the  conditions  of  this  experiment. 
Values  of  rs  of  14.0  sec/cm  or  greater  make  the  latent  heat  term  of 
equation  (77)  negligible. 
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TO 

Results  and  Discussions 

The  following  equation  can  he  derived  (Gebhart  (l96l))  from  the 
Folhausen  similarity  solution  for  the  two-dimensional  boundary  layer  of 
a  flat  plate. 

ra^  =  1.92  (L'/u)l/2  Pr2/3  (78) 

L'  is  the  width  of  the  plate,  u  is  wind  velocity  and  Pr  is  the  Prandtl 
number.  The  Prandtl  number  is  the  kinematic  viscosity  of  the  fluid  (in 
this  case,  air)  divided  by  its  thermal  diffusivity .  Figure  (24)  shows 
measurements  of  ra^  plotted  against  (L' /u)1/2(Pr)2/3  for  a  horizontal 
leaf.  Six  to  ten  measurements  at  different  radiation  loads  were  made  at 
a  given  wind  velocity  for  a  particular  leaf.  No  significant  correlation 
was  found  between  the  leaf  to  air  temperature  differences  and  ra  at  even 
the  lowest  wind  velocity  (4l  cm/sec).  Because  of  the  difficulty  in  main¬ 
taining  a  constant  air  temperature  at  the  lowest  wind  velocities,  the 
method  was  not  precise  enough  to  study  natural  or  mixed  convection.  Thus, 
means  and  standard  deviations  were  calculated  from  the  six  to  ten 
measurements  at  each  wind  velocity  for  a  particular  leaf.  A  least  squares 
method  was  used  to  find  the  best  linear  fit  to  the  data  using  all  the 
measurements.  The  line  was  forced  through  the  origin.  The  sum  of 
squares  of  the  deviations  from  this  line  was  used  to  calculate  a  sample 
standard  deviation  of  the  slope  which  was  in  turn  used  with  the  ’ t ’ 
distribution  (Snedecor  (1956))  to  establish  a  confidence  interval  for 
the  slope  at  the  95%  level.  The  slope  and  confidence  interval  for  the 
data  in  Figure  (24)  is  1.51  ±  0.032  sec1//2/cm.  This  is  significantly 


lower  than  the  slope  of  the  theoretical  line.  Since  a  leaf  has  leaf  hairs 
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.  Boundary  layer  resistance  for  heat  ra^  as  a  function  of 

(L'/u)1'2(Pr)  for  a  leaf  in  the  horizontal  position.  L'  (cm) 
is  the  leaf  width,  u  (cm/sec)  is  the  wind  velocity  and  Pr  is 
the  Prandtl  number. 


ridges  undulations  etc. ,  it  is  not  surprising  that  a  leaf  has  a  boundary 
layer  resistance  lover  than  that  of  a  flat  plate. 

Figure  (25)  shews  similar  measurements  for  the  leaf  in  a  natural 
or  more  upright  position.  This  slope  with  95$  confidence  interval  is 

*i/p 

1.72  ±  0.01*0  sec  '  /cm.  These  values  are  significantly  higher  (the 
confidence  intervals  do  not  overlap)  than  corresponding  values  for  the 
leaf  in  the  horizontal  position.  Die  leaf  in  the  upright  position  acts 
as  a  bluff  body  in  the  air  stream.  Separation  occurs  as  the  air  moves 
around  the  leaf.  The  movement  of  air  is  retarded  close  to  +he  underside 
of  the  leaf.  Thus  heat  transfer  from  the  underside  is  less  than  that  for 
the  side  facing  the  air  stream.  This  effect  of  separation  on  heat  trans¬ 
fer  per  unit  area  has  been  verified  experimentally  with  cylinders  (Gebhart 
(1961)).  This  effect  would  be  especially  noticeable  for  corn  leaves  in 
this  upright  position  because  the  edges  of  the  leaf  tend  to  fold  around 
the  mid  rib  creating  a  half  cylinder  effect.  This  folding  could  lead  to 
appreciably  less  heat  transfer  per  unit  area  than  a  leaf  at  zero  angle 
of  attack. 

The  effect  of  free  stream  turbulence  on  the  boundary  layer  resis¬ 
tance  is  shown  in  Figure  ( 26) .  With  no  leaves  in  the  wind  tunnel  the 
heated  thermocouple  annemometer  detected  less  than  1.0#  turbulent 
intensity.  When  two  leaves  were  placed  upwind  of  the  annemometer  tur¬ 
bulent  intensities  ranged  from  a  maximum  of  15-5 #  at  86  cm/sec  to  a  mini¬ 
mum  of  5.^#  at  330  cm/sec.  The  slope  of  the  best  linear  fit  to  the 
measurements  in  Figure  (26)  with  its  95 #  confidence,  interval  is  1.^7  ± 
0.0^3  sec  1/1,2 /era.  Thus,  the  turbulence  generated  by  two  leaves  upwind 
caused  a  15#  decrease  in  boundary  layer  resistance.  Ibis  result  is 
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Figure  25.  Boundary  layer  resistance  for  heat  (ra^)  as  a  function  of 
(L’ /u)1/2(Pr)  for  a  leaf  in  an  uprigrt  position.  L'  (cm) 
is  the  leaf  width,  u( cm/sec)  is  the  wind  velocity  and  Pr 
is  the  Prandtl  number. 


Figure  26.  Boundary  layer  resistance  for  heat  (ra^  as  a  function  of 
(L'/u)1/2(Pr)  for  a  leaf  in  an  uprigit  position  with  two 
leaves  upwind.  L'  (cm)  is  the  leaf  width,  u  (cm/sec)  is  the 
wind  velocity  and  Pr  is  the  Prandtl  number. 
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comparable  to  those  of  Maisel  and  Sherwood  (1950 ).  Working  with  spheres 
at  Reynolds  numbers  similar  to  this  study  (5,000  to  20,000)  they  found 
a  25 %  decrease  in  boundary  layer  resistance  when  turbulent  intensities 
were  increased  from  h%  to  25$. 

The  difficulty  remains  as  to  how  to  extrapolate  wind  tunnel  measure¬ 
ments  to  the  field.  Obviously,  the  structure  of  turbulence  in  the  field 
is  quite  different  than  the  free  stream  turbulence  generated  by  leaves 
in  the  wind  tunnel.  However,  this  wind  tunnel  study  serves  as  a  useful 
comparison  with  field  measurements.  Data  from  Kanemasu  et  al.  (1969 ) 
obtained  from  measurements  made  under  field  conditions  on  attached 
bean  leaves  are  shown  in  Figure  (27).  Hie  measurements  are  plotted 
against  (L/u)1^.  The  increase  of  turbulent  exchange  from  the  leaf 
surface  in  the  field  should  lead  to  a  suppression  of  the  Prandtl  number 
effect.  Thus  functions  of  (L.A’)l/2  rather  than  (L' /u)1//2(Pr)2^  will 
serve  for  field  conditions.  This  also  means  that  ra^,  ray.  and  rac  are 
assumed  equal.  The  slope  with  its  95%  confidence  limit  of  the  least 
squares  fit  to  the  data  of  Kanemasu  et_  al.  is  0.60  ±  0.31  sec^^/cm  for 
a  5  cm  leaf  width  and  0.1+3  ±  0.19  sec^/^/cm  for  a  10  cm  leaf  width.  The 
data  in  Figure  (27)  are  plotted  using  the  5  cm  width.  Also  shown  in 
Figure  (27)  are  the  best  fit  to  the  wind  tunnel  data  with  and  without 
leaves  upwind.  The  slope  is  reduced  by  approximately  15$  by  turbulence 
generated  by  two  leaves  in  a  wind  tunnel.  The  field  measurements  repre¬ 
sent  a  further  reduction  of  about  56$. 

Because  of  lack  of  knowledge  of  the  nature  of  turbulence  immediately 
above  and  within  a  plant  community,  one  can  only  speculate  as  to  causes 
of  the  relatively  low  boundary  layer  resistances  in  the  field.  However, 
very  high  turbulent  intensities  have  been  reported  under  field  conditions. 


(Allen  (1968)).  Also  movement  of  leaves  and  unsteady  wind  conditions 
would  result  in  transient  boundary  layers  and  disrupted  flow  around  the 
leaf.  Also,  speculations  have  been  made  (Hunt  at  al.  1966)  on  the 
existence  of  turbulent  boundary  layers  on  leaves  in  the  field.  How¬ 
ever,  turbulent  boundary  layers  in  the  classical  fluid  mechanics  used 
of  this  word  exist  at  Reynolds  numbers  much  higher  than  those  encountered 
in  this  study.  Figure  (27  )  shows  the  theoretical  curve  for  a  turbulent 
boundary  layer  for  a  flat  plate  (Gebhart  (l96l)).  The  turbulent  r  values 
become  abnormally  large  as  (L'/u)^^  increases  because  the  theory  does  not 
apply  at  low  wind  speeds.  A  review  article  by  Kestin  (1966)  makes  clear 
that  experimental  measurements  verify  the  theory  at  high  Reynolds  numbers. 
An  extrapolation  of  this  theoretical  and  measured  curve  to  lower  Reynolds 
numbers  intersects  the  curve  describing  the  laminar  boundary  layer  at 
a  Reynolds  number  of  about  50,000.  (The  highest  Reynolds  number  used  in 
this  study  was  approximately  20,000.)  However,  it  is  clear  that 
free  stream  turbulence  does  enhance  heat  transfer  at  relatively  low 
Reynolds  numbers.  Considering  the  high  turbulent  intensities  and  unsteady 
state  conditions  encountered  in  the  field,  measurements  of  Kanemasu  et  al. 
(1969)  represent  a  reasonable  extrapolation  of  the  wind  tunnel  measure¬ 
ments  reported  in  this  study. 
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Conclusions 

4 

Measurements  of  the  boundary  layer  resistance  for  heat  (ra^)  made 
in  a  wind  tunnel  can  be  described  by  the  following  equation; 

rah  =  1.72  (L'/u)l/2(Pr)2/3  (79) 

Free  stream  turbulence  generated  by  two  leaves  upwind  of  a  third  leaf 
caused  a  reduction  of  ra^  on  the  third  leaf  of  approximately  15$. 

Thus,  ra^  under  these  conditions  can  be  described  by 

re ^  =  1.51  (L'/u)l/2(Pr)2/3  (80) 

Turbulent  intensities  measured  downwind  from  the  two  leaves  ranged  from 
15.5#  at  86  cm/sec  to  5-4#  at  330  cm/sec.  Under  field  conditions, 
measurements  by  Kanemasu  et  al.  (1969 )  suggest  the  following  function, 

ra  =  0.60  (L'/u)l/2  .  (81) 

This  equation  represents  a  further  reduction  in  ra  values  of  56#. 
Unsteady  state  and  highly  turbulent  conditions  in  the  field  suggest  that 
the  function  derived  from  the  measurements  of  Kanemasu  et  al.  represents 
a  reasonable  extrapolation  of  the  wind  tunnel  measurements  reported  here. 
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Chapter  IV 


SIMULATING  THE  CROP  RESPONSE 

The  net  photosynthesis  rate  of  an  individual  leaf  was  expressed  in 
Chapter  II  as  a  function  of  the  visible  radiation  load,  the  leaf  temper¬ 
ature  and  the  C0£  concentration  at  the  leaf  surface.  Presumably,  if  these 
variables  can  be  estimated  for  all  the  leaves  in  a  plant  community,  the 
individual  leaf  photosynthesis  rates  can  be  integrated  over  the  total 
leaf  area  to  obtain  the  photosynthesis  rate  of  the  crop.  A  method  of 
calculating  the  distribution  of  visible  radiation  was  described  in 
Chapter  I.  Discounting  the  effects  of  movement  of  leaves,  this  distri¬ 
bution  of  visible  radiation  is  independent  of  the  air  stream  above  and 
within  the  canopy  of  plants.  The  leaf  temperature  and  CO  concentration 
at  the  leaf  surface  are  dependent  on  the  air  stream.  The  problem  becomes 
one  of  describing  in  mathematical  terms  the  turbulent  transport  of  heat, 
momentum  and  mass  from  some  reference  height  above  the  crop  surface  to 
the  individual  leaf  surfaces.  The  problem  actually  has  two  boundaries 
since  conditions  both  at  this  reference  height  and  at  the  soil 
surface  must  be  known.  This  system  of  equations  which  describe  the 
turbulent  transport  between  the  leaf  and  the  environment  must  of  necess¬ 
ity  include  source  and  sink  terms  caused  by  the  amount  of  CO,-.,  sensible 
and  latent  heat  given  off  or  taken  up  by  the  leaves  themselves.  That  is, 
the  leaf  temperatures  and  CO^  concentrations  are  dependent  both  on  the 
exchange  properties  of  the  air  stream  and  the  leaf  characteristics 
themselves.  This  interaction  can  be  studied  with  such  a  system  of 
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equations .  The  method  of  formulating  turbulent  transport  and  the 
solution  of  this  system  of  equations  which  include  the  calculations  of 
the  rates  of  net  photosynthesis  as  well  as  sensible  and  latent  heat 
exchange  will  now  be  considered  in  detail. 
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Theoretical  Consideration 

Above  the  Vegetation 

Many  of  the  theoretical  treatments  to  date  have  placed  the  upper 
boundary  at  or  near  the  top  of  the  crop  (Philip  (196U),  Waggoner  and 
Riefsnider  (1965).  Ccwen  (1968)).  However,  the  top  of  the  crop  is 
difficult  to  define  in  terms  of  height.  Also  the  variation  of  wind 
velocity,  air  temperature,  etc.  with  height  is  at  a  maximum  in  this 
region.  Thus  a  better  boundary  to  the  problem  is  at  an  arbitrary  height 
say  two  to  four  meters  above  the  crop,  yet  still  within  the  crop  boundary 
layer,  where  the  variables  such  as  air  temperature  and  COg  concentration 
are  approximately  constant  and  where  the  shearing  stress  and  fluxes  from 
the  crop  are  constant  with  height.  In  this  study,  four  meters  above  the 
crop  was  chosen  as  the  reference  height. 

The  exchange  processes  for  mass,  heat  and  momentum  are  described  by 
one-dimensional  continuity  equations  both  above  and  within  the  plant 
community.  For  momentum 

p-=  0  (80) 

3z 

where  t  =  pKM  (8l) 

t  is  the  shearing  stress,  p  is  the  air  density,  is  the  eddy  diffusiv- 
ity  for  momentum,  u  is  the  wind  velocity  and  z  is  the  distance  in  the 
vertical  direction.  From  considerations  of  the  Reynolds  stresses 
and  the  Prandlt  mixing  length  concept  (Sutton  (1955)),  the  neutral  wind 
profile  can  be  derived  as ; 
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u  =  ^  ln(i^)  (83) 

z0 

where  u*  is  the  friction  velocity  and  is  expressed  as; 

u*  =  A/p 

Other  undefined  variables  are  k,  the  von  Karmen's  constant  (taken  to  be 
0.4  in  this  treatment),  D  the  zero  plane  displacement  and  zD  the  roughness 
length.  A  crop  forces  the  reference  plane  upwards  from  the  ground  surface 
to  D.  z0  is  the  distance  above  D  where  the  wind  profile  extrapolates 
to  zero.  In  this  treatment,  zQ  and  D  are  treated  as  crop  constants 
and  are  estimated  from  measurements  of  z0  and  D  in  com  as  functions  of 
the  crop  height  which  were  presented  in  Lemon  (1965). 

For  diabatic  conditions,  the  theory  developed  by  Swinbank  (1964) 
and  the  Keyps  formulation  (Ellison  (1957).  Panofsky  (1961),  Sellers  (1962)) 
are  most  commonly  used.  In  differential  form  the  Swinbank  and  Keyps 
equations  which  relate  wind  velocity  to  height  are  given  respectively  as 


*u=  f  (1  -  exp(  -  :  D)  ir1 

a  z  kz  L 

(84) 

and 

sk  +  x  — r~)s^  ' 1 

(85) 

where 

„  _  k(z  -  D)  9u 

k  u*  9z 

(86) 

and 

u#3T  P  CP 

L  =  -  r 

(8?) 

k  g  H 

H  is  the  sensible  heat  flux  from  the  crop  surface  and  is  considered 
positive  in  the  upward  direction.  (This  sign  convention  is  also  used  for 


net  photosynthesis  and  latent  heat  flux).  L  is  the  Monin-Obukhov  length, 

X1  is  a  constant  and  is  the  eddy  diffusivity  for  heat  and  is  defined 
as 

KH  =  H/fcCp|i)  (88) 

where  T  is  the  air  temperature.  Both  of  the  above  theories  have  been 
modified  to  describe  the  exchange  of  heat  and  mass  across  the  crop  boun¬ 
dary  layer.  For  the  purposes  of  this  study,  both  were  in  general  agreement. 
Hie  Swinbank  approach  is  incorporated  into  this  model  because  it  is 
basically  simpler  and  involves  much  less  numerical  computation.  Thus, 
only  this  approach  will  be  described  here. 

By  integrating  (8k)  one  obtains 


(exp(zo)  -  l)  ) 
L 


Combining  (8l)  with  ( 8U )  one  obtains 


(89) 


KM  =  u*  k  L  (1  -  exp(-z,/L) }  (90) 

Swinbank  (1968)  presented  evidence  that  the  eddy  diffusivity  for  heat 
Kjj  and  momentum  (Km)  are  related  by 

jT  “  2.7(z/L)0-21*  (91) 

km 

for  unstable  conditions.  There  are  two  difficulties  with  this  expression. 
It  obscures  the  fact  that  the  ratio  K^/Kj^  approaches  a  constant  value 
as  (z/L)  increases  (That  is  as  the  free  convection  region  is  approached 
(Ellison  (1957))*  Also  at  low  values  of  |  z/L )  (near  neutral  conditions) 


this  ratio  falls  as  low  as  1.2  or  lower.  Data  presented  by  Panofsky 
(1965)»  (Figure  (28))  indicates  that  this  ratio  is  much  higher  at  near 
neutral  conditions.  Therefore  K^/K^  was  defined  as; 


%/%  =  -l.U  exp(l.5(z/L))  +  3.0 


Using  equation  f'2),  %/%  approaches  3.0  as  |  z/L|  approaches  infinity 
and  at  neutral  conditions  (z/L)  becomes  1.6.  Figure  (28)  shows  both 
equations  (91)  and  (92)  along  with  data  plotted  by  Panofsky  (1965). 
These  data  were  taken  from  studies  by  Swinbank  (1961*). 

Combining  equations  (88)  and  (90)  and  (92)  one  obtains 


3T  H  3 


3z  pCpU*^(-l.l*  exp(l.5(z/L) )  +  3)  (1  -  exp(-z/L))  3z  (93 


Integrating,  one  obtains 


Tch  Th  ~ 


pCpU*L 


(-1.1*  exp(-i. 5y )  +  3)  (l  -  exp(y)) 


=  Th  - 


pCpU*L 


where  y  =  (-z/L), 

T  ,  is  the  air  temperature  at  a  height  (ch)  immediately  above  the  crop 

Cli 

and  Tjj  is  the  air  temperature  at  the  reference  hei$it  (h)  four  meters 
above  the  crop.  FY  is  the  integral  in  equation  (9I*).  Similar  expressions 
for  water  vapor  pressure  (e)  and  CO2  concentration  (c)  are  respectively, 

ech  “  eh  pCPu*kL  (96) 


and  C.v,  =  C 


86 


LE  and  are  the  latent  heat  flux  and  the  total  net  photosynthesis  of 
the  crop.  Equations  (96)  and  (97)  are  derived  by  assuming  the  eddy 
diffusivities  for  heat  and  mass  transfer  are  identical.  Wind  velocity, 
the  eddy  diffusivity  for  momentum  and  the  eddy  diffusivity  for  sensible 
heat  at  the  top  of  the  crop  are  expressed  respectively  as. 


uch  =  ”h 


—  i  j 
-  In  < 


exp(^_ D.)  -  1  ") 

exp (-'■  -)  -  1  ) 


(98) 


KMc  =  u*  k  L  {1  -  exp(^—  £) } 

L 

and  =  K^c(— l.U  exp(l.5(ch  -  D)/L)  +  3) 


(99) 

(100) 


Equations  (9M  -  (100)  assume  that  u*  and  L  are  known.  They  can 
be  determined  if  z0,  D  and  H  are  known.  The  two  equations  with  u*  and 
L  as  unknowns  are. 


u*  =  k  l^/ln 


expfo-  ~  -  1 

L 

expLial  -  1  ' 

L 


(101) 


and 


u*^pCp 

L  =  - 

k  g  H 


(102) 


They  can  be  solved  numerically  by  the  method  of  variable  position 
described  in  Conte  (1965).  Therefore,  the  above  theory  allows  us  to 
estimate  the  needed  variables  at  or  near  the  top  of  the  crop  from  measure¬ 
ments  at  some  reference  height  above  the  crop.  Note  that  the  crop  para¬ 
meters  D  and  z0  and  the  values  of  H,  and  LE  must  be  known.  This 
interrelationship  between  the  changes  in  the  air  properties  and  the  sources 
and  sinks  will  be  explained  in  detail  below. 
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Exchange  within  the  Crop 

In  the  crop,  the  equations  of  continuity  for  momentum  and  mass  and 
sensible  heat  transfer  are; 


3t  = 
dz 

1  pu2  C,  f 

2  d 

(103) 

and 

3£  = 
3z 

Q"f 

(IOU) 

Where  Q  represents  the  vertical  flux  of  mass  or  sensible  heat  and  Q" 
represents  the  flux  of  the  same  quantities  from  a  unit  area  of  leaf 
surface.  Cd  is  a  drag  coefficient  and  f  is  the  leaf  area  density.  Follow¬ 
ing  Sutton  (1955),  the  shearing  stress  can  be  expressed  as 

t  =  pu'v'  (105) 

where  u'  and  v'  are  deviations  from  the  mean  horizontal  and  vertical 
wind  components.  Assuming  u'  is  equal  to  v'  and  that 

u'  =  £—  (106) 

dz 

x  can  be  expressed  as; 

x  =  £2  —  (107) 

dz 

Therefore,  equation  ( 10 3 )  can  be  expanded  to 


3_ 

9z 


<%  !x>  -  1/2  0Cd  u*f 


(108) 


h  (p£2(!l)2)  -  1/2pc*  u2f 


(109) 


Equations  (108)  or  (109)  are  usually  solved  by  assuming  £  or  and  f 
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are  constants  or  linear  functions  with  height  (Cionco  (1965),  Lemon  et  al. 
(1963)).  Assumptions  about  f  can  be  avoided  by  changing  the  independent 
variable  from  z  to  F,  the  cumulative  leaf  area  index  measured  from  the 
top  of  the  crop  downward.  (Perrier  and  Hallaire  (1968)).  The  transforma¬ 
tion  is  given  by, 

M  =  -f 

3z  (110) 

Equation  (103)  can  be  solved  by  assuming  the  turbulent  intensity  (b) 
is  constant  with  L  or  z.  The  turbulent  intensity  is  expressed  as 


u 


(111) 


Combining  (103),  (105)  and  (ill)  results  in 

B2^  =  1/2  Cdu2  f  (112) 

3z 

Then  u  =  exp(  -  F) 

2B2 

or  u  =  Uch  exP(-f'°F )  (H3) 

where  y0  is  a  constant.  The  difficulty  with  this  equation  is  that  the 
region  near  the  top  of  the  crop  has  very  little  leaf  area  but  is  where 
the  wind  velocity  decreases  a  relatively  large  amount.  To  apply  equation 
(113)  to  this  region  would  produce  unreasonably  large  values  of  yQ.  In 
this  same  region,  the  exponential  wind  equation  (89)  predicts  the  wind 
profile  reasonably  well.  Thus  a  height  (h^)  with  a  cumulative  leaf 
area  index  of  F^  below  the  top  of  the  crop  but  above  D  +  z0 ,  was  chosen 
where  the  wind  profile  begins  to  deviate  from  the  exponential  law.  Then 

u  =  Uhl  exp(-yo(F  -  F^)  (llU) 
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was  used  to  predict  the  wind  profile  further  into  the  crop  where 


«hl  =  ^ 


-T 


exp(*— ^)  -  1 
_ L _ 

exp  D  -  1 

L 


(115) 


The  constants  y0  and  h-^  were  determined  empirically  from  wind  profile 


data. 


Hie  ed<3y  diffusivity  in  the  crop  was  assumed  to  be  proportional 
to  the  wind  velocity.  Cowan  (1968)  made  this  same  assumption.  Hie  pro¬ 


portionality  factor  (6)  was  determined  by 


(116) 


Hie  ratio  of  the  eddy  diffusivities  for  heat  and  momentum  in  the  crop  was 
assumed  constant.  Wright  and  Brown  (1967)  have  shown  that  this  is  a 
reasonably  good  assumption.  Thus  in  the  crop 


(117) 


KHc 

h  *  <ir° 

Me 


Kjj  =  (-Ha)  6u 

^c 


(118) 


Again  the  eddy  diffusivities  for  mass  transfer  are  assumed  equal  to  K^. 
Sensible  Heat  and  Mass  Exchange 

The  transfer  of  sensible  heat  and  mass  transfer  within  the  crop  is 
described  by  equation  (104).  For  sensible  heat,  equation  (104)  can  be 


expressed  an ; 


(119) 
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where  is  the  leaf  temperature  and  T  is  the  air  temperature  of  the 
bulk  air  stream  surrounding  the  leaf.  In  actual  fact,  at  a  given  level 
in  the  crop,  leaves  will  be  at  a  range  of  temperatures  governed  mainly 
by  their  respective  radiation  loads.  In  Chapter  I  a  method  was  presented 
which  divided  the  leaf  area  of  a  given  leaf  area  increment,  into  classes 
of  visible  radiation  loads.  Ihis  same  procedure  can  be  used  to  divide 
the  same  leaf  area  into  classes  of  solar  radiation  loads.  This  procedure 
will  be  discussed  in  more  detail  below.  Since  twenty  classes  were  used, 
equation  (119)  can  be  written  as 


20 

A  (%j  -  T*> 

-3T=0CP  ^ - 


rQ. 


FR^ 


(120) 


FR^  represents  the  frequency  distribution  of  the  leaf  area  classes  at 
the  leaf  area  increment  i.  Since 


H  =  +f 


pCp 


3T 

8L 


equation  (120)  can  be  expressed  as; 


3_ 

3L 


,  9Ti  X 

(fiC3T)  = 


Ti> 


ra. 


(121) 


(122) 


Replacing  the  differentials  with  finite  differences,  one  can  write 


-  fiKi 


(Ti-1  -  Ti> 


+  fi+l  Ki+1 


(Ti  -  W  ^TL-m  -  Ti) 


=  S 


±LL 


FRij  (123) 


where  S  is  the  leaf  area  per  increment.  Figure  (29)  represents  a  crop 
divided  into  n  increments.  Note  that  is  the  air  temperature  at  the 


reference 

height 


Figure  29.  A  representation  of  the  crop  divided  into  n  leaf  area  incre¬ 
ments.  Air  temperature  and  leaf  area  density  are  shown  as 
functions  of  height  (z). 
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center  of  inc  .‘ment  i.  On  the  other  hand,  f^  and  are  estimated  at 
the  dividing  line  between  increment  i-1  and  increment  i. 

The  subscripted  K’s  refer  to  the  eddy  diffusivity  for  mass  and  heat 
transfer  and  can  be  evaluated  from  equation  (118).  The  values  of  r_  can 
be  evaluated  from  (llM  and  the  relation 


ra  =  0.6  (L'/u)l/2 


(12  It) 


determined  from  measurements  of  Kanemasu  (1968)  which  were  presented  in 
Chapter  III.  L'  is  the  leaf  width.  Equation  (123)  can  also  be  written  as 

-<fiKi>Ti-i  +  <Vi  ♦  ♦  ->Ti  -  «WWTn-i  ■ 


- —  E  T  FR. 

r*i  h}  'i 


(125) 


The  corresponding  equation  at  the  first  increment  just  below  the  top  of 
the  crop  is ; 


-  TX  Tx  -  T2  E(Tl  -  T-.) 

-  )  +  Klfl(  - )  =  S  _ ij _ L  FR  (126) 

Az  S  rn.  J 


.  Tch  -  Tl  % 

-Kch(  .  )  +  Klfl( 


Az  is  chosen  so  that  T^  is  the  air  temperature  at  the  middle  of  the  first 
leaf  area  increment.  That  is,  Az  is  twice  the  distance  from  the  middle 
of  the  first  increment  to  the  top  of  the  crop.  Fortunately  the  change 
of  temperature  with  height  is  quite  linear  in  this  region.  Thus,  rela¬ 
tively  large  values  of  z  do  not  introduce  appreciable  error.  At  the 
increment  of  leaf  area  next  to  the  soil  surface,  equation  (123)  becomes 


-Vn 


(Tn-1  “  ^n) 


-  FHS  =  -  E  T. 

Pan 


Lnj 


(127) 


where  FHS  represents  the  flux  of  sensible  heat  from  the  soil.  Thus, 
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assuming  leaf  temperatures,  FHS  and  are  known,  equations  ( 125 )  (126) 
and  (127)  represent  a  system  of  n  linear  equations  in  n  unknowns. 

This  system  of  equations  can  best  be  treated  in  matrix  notation 
expressed  as 

(AM  +  BM)  (T)  =  (FH)  (127) 

or  (T)  =  (FH)  (AM  +  BM)-1 


where 


■4'f 

IT 


9b 


!"l,  ,  '  rau  ‘  s2/r»i  *  K„>,  *  •  S/  z 


•  v  -  #  # 

FH  = 


<4*.*  -r**  *•♦  -•  ► 


ch  ch 


STr«,  *  FH?1  •  S2/ 


L2J  *  ™2J  ‘  s  'ra2 


rrL  •  FRjjj  •  s3/^  +  fhs 


nj 


and 


T  = 


A  similar  equation  for  water  vapor  is  E 


E  =  (FLE)  (AM  +  EM) 


-1 


f  S2  EFR. ,  / (rai  +  rs-y) 


EM  = 


i j ' v  al 
0 


S2  ZFR2j/(ra2  +  rs2j) 


0 
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S2  +  rs13)  *  S/4z 

S2  ^e(T^)/(ra2  t.rs^) 

S2  Ze(T^)/(raii  +  rsnj)  +  FLHS 
where  FLHS  is  the  flux  of  latent  heat  from  the  soil 


where  ei  is  the  water  vapor  pressure  of  the  air  stream  at  the  i*11  incre¬ 
ment.  Die  values  of  the  resistance  rs^j  are  derived  from  equation  (48) 
expressed  in  subscripted  form  as; 


is  the  light  flux  density  at  class  j  of  increment  i.  Since  the  leaf 
area  classes  for  visible  radiation  loads  were  calculated  simultaneously 
with  the  leaf  area  classes  for  solar  radiation  loads,  it  was  not  diffi¬ 
cult  to  assign  a  leaf  resistance  value  to  each  class  of  solar  radiation 
loads . 

The  equation  describing  C02  exchange  is  solved  in  a  less  direct 
manner.  Net  photosynthesis  (N)  per  unit  leaf  area  is  related  to  CL,  the 
C02  concentration  at  the  leaf  surface  by  equation  (57).  (This  equation 
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will  be  written  here  simply  as  N  =  fn(CL)).  The  continuity  equation  for 


C02  exchange  can  be  written  as 

-  ---ijr." - - 

■  Tp =  fn<CL) 


(131) 


where  q  is  the  vertical  flux  of  COg.  Once  again  there  is  a  range  of 
N's  and  CL's  in  a  given  leaf  area  increment  which  depend  on  the  distri¬ 
bution  of  visible  radiation  loads  and  are  represented  by  N.  .  and  CL.  , . 

ij  ij 

The  leaf  surface  concentrations  of  COg  can  be  expressed  as 


CL  =  q  -  raiNit3 


Thus  equation  (131)  can  now  be  written  sis 


=  Zfh(q  -  rsujN^ )  Fqj 


(132) 


(133) 


In  matric  notation  this  reduces  to 


(AM) (C)  =  FC 


where 


S2  Slhtq  -  ra^jFEy  +  S/ 


(134) 


(0  = 


S2  2fn(  q  -  r^qjjFRnj  + 


.  V,. 
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FS  and  are  the  flux  of  CO^  from  the  soil  and  the  CC>2  concentration 
at  height  ch  respectively.  In  this  situation  it  was  impossible  to 
separate  the  air  CO^  concentrations  into  a  separate  vector.  Thus  equation 
(13^)  was  solved  iteratively  for  (C).  For  the  initial  guess,  all  C's 
were  set  equal  to  the  concentration  of  CO2  at  the  reference  height. 

Using  matrix  theory  and  in  the  case  of  CO^  an  additional  iteration 
sequence,  the  vectors  (T),  (E),  and  (C)  are  evaluated.  Once  they  are 
known  the  sources  and  sinks  of  sensible  (SO^)  and  latent  (SOg)  heat  and 
CO  (soc)  axe  calculated  for  each  leaf  area  increment  from; 

S0Ri  =  pCp  KT-  -  Ti)  FRij  S/rai 

S°Ei  =  PCp  ECeCTj^j)  -  ej)  FR^  S/(rai  +  rs^ ) 
and  SO^  =  ZfnCC^  -  rail^ . )  FR^  S 

The  vertical  flux  of  sensible  heat  from  the  crop  surface  is  obtained  by 
summing  the  sources  and  sinks  over  all  the  leaf  area  increments  and 
adding  the  flux  from  the  soil  surface.  A  similar  procedure  is  used  for 
latent  heat  and  COg. 

Leaf  Temperatures 

The  leaf  temperatures  are  obtained  by  solving  the  energy  balance 
equation  for  individual  leaves.  This  is  expressed  as 


Rn.  =  H.  +  LF. 


(135) 


where 


Rni^  =  1%+  +  THi+1+  -  2aTL.j  +  Aj  •  DS  *0.3/sin(lS). 


Hij  =pCP  (TLU  -  Ti)/rai 
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PC  e(T  )  -  ei 

LE  =  -L  (—hj - i-) 

Y  rs^  +  r&i 


(136) 


Aj  is  a  fraction  of  twenty  minus  0.025.  That  is,  is  one  of  0.05  - 
0.025,  0.10  -  0.025,  0.15  -  0.025,  ...  1.0  -  0.025.  The  number  0.025 
is  subtracted  to  put  the  radiation  flux  density  in  the  middle  of  its  res¬ 
pective  class.  These  flux  densities  Eire  expressed  as  fractions  of  DS  •  0.7/ 
sin (IS)  which  represents  approximately  the  largest  solar  radiation  load. 

DS  is  the  amount  of  direct  solar  radiation  measured  on  the  horizontal 
plane  above  the  crop.  TH^t  and  TH^+^+  are  the  amounts  of  theriml  radia¬ 
tion  moving  up  and  down  in  the  crop  as  discussed  in  Chapter  I.  An 
anomaly  arises  here  because  the  leaf  temperatures  must  be  known  before 
the  thermal  radiation  values  can  be  calculated.  This  difficulty  is  dis¬ 
cussed  below  in  connection  with  successive  approximations.  Equation 
(135)  represents  20  non-linear  equations,  for  each  increment  i.  To 
save  computation,  only  every  sixth  equation  was  solved  for  Ty  using 

ij 

once  again  the  method  of  variable  position.  3he  leaf  temperatures  for 
in-between  classes  were  found  by  linear  interpolation.  Also  the  follow¬ 
ing  equation  was  used  to  adjust  the  leaf  temperatures. 


9Rn  _  9H 
3F  9F 


9Rn  _  n 
3F  p°p 


^  -  *i)raiJ  pC  (eUT^)  -  ei)FR.J 


rai  +  rsij 


Since  net  radiation  as  a  function  of  F  is  calculated  in  the  model  the 
above  equation  can  be  solved  iteratively  for  i)>.  Then  was  used  to 
adjust  the  leaf  temperatures  so  that  the  SRjj/SL  was  E'qual  to  the  sources 


of  sensible  and  latent  heat  for  a  given  increment. 
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The  Soil  Surface 

The  soil  surface  represents  the  second  boundary  of  this  two  boundary 
mathematical  problem.  In  this  treatment,  the  soil  surface  is  defined 
in  terms  of  a  heat  flux  into  the  soil  (SHF)  and  a  soil  moisture  tension  (SM). 
The  first  term  is  measured  directly  using  soil  heat  flux  plates  near  the 
surface.  The  second  is  estimated  from  soil  moisture  tension  measure¬ 
ments  made  at  various  depths  in  the  soil.  Mass  and  heat  exchange  between 
the  soil  surface  and  the  air  stream  above  the  soil  is  described  using  the 
theory  of  Owen  and  Thompson  (1963).  This  approach  has  been  used  exten¬ 
sively  by  Chamberlain  ( 19 68)  who  studied  the  movement  of  gases  to  and 
from  a  variety  of  rough  surfaces . 

The  resistance  to  diffusion  from  the  surface  to  the  air  stream  is 


(137) 

(138) 


where  Hs  and  LES  are  fluxes  of  sensible  and  latent  heat  per  unit  area 
of  soil  respectively.  TS  is  the  soil  temperature,  eg  is  the  vapor  pressure 
at  the  soil  surface,  ezh  and  T^  is  the  air  water  vapor  pressure  and 
air  temperature  at  a  reference  height  (zh)  above  the  soil  surface.  This 
was  arbitrarily  taken  as  the  middle  of  the  n  leaf  area  increment.  In 
these  equations,  the  friation  velocity  u*  is  defined  as 


uzh  k 
ln(z/zQ) 


(139) 


where  is  the  velocity  of  the  air  stream  at  the  reference  height  and 
zQ  is  a  roughness  length  which  is  a  characteristic  of  the  soil  surface. 

■*  -  .-Ike  .jres ietenoss-  were  determined*  from  the*:follawirig*eqtC&,'tiOn‘by  •  "  ' 

Owen  and  Thompson  (1963), 

rzH  =  0.52(  ^  (Pr)°-8  (1U0) 

U#  V 


where  v  is  the  kinematic  viscosity  of  air  and  zh  is  the  reference  height. 
For  rzE,  the  Prandtl  number  (Pr)  is  replaced  by  the  Schmidt  number,  foe 
various  numerical  constants  in  equation  ( lU 3 )  were  determined  by  measure¬ 
ments  of  Chamberlain  (1966). 

foe  energy  balance  for  the  soil  surface  can  now  be  expressed  as 


'  ns 


pCP  ,es  -  Ts  Tzh 

+  PV1*  - )+  SHF 


raH 


Uhl) 


The  vapor  pressure  at  the  soil  surface  is  derived  from 


es 


e(Ts)  exp( — -^— ) 
R  T 


v  s 


(1U2) 


where  e(Ts)  is  the  saturated  vapor  pressure  at  temperature  Tg  and  Ry  is 

the  gas  constant  for  water  vapor.  R_  is  defined  as; 

“s 


Rng  =  VIS(1  -  ALB1)  +  NIR(l  -  ALBg)  +  TH+  (n)  -  (TS)1*  (lU3) 


where  VIS  and  NIR  are  the  total  amounts  of  visible  and  near  infrared 
radiation  respectively  which  penetrates  to  the  bottom  of  the  canopy; 

ALBj  and  ALB2  are  the  soil  albedos  for  visible  and  near  infrared  radiation 
and  TH+(n)  is  the  thermal  radiation  moving  downward  from  the  n^*1 


increment  of  leaf  area.  Once  again,  a  non-linear  equation  must  be  solved. 
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Again  we  solve  equation  (li+5)  for  Tg  using  the  method  of  variable  position. 
Once  T  is  known,  the  fluxes  of  sensible  and  latent  heat  from  the  soil 
surface  can  be  calculated  using  equations  (137)  and  (138)  respectively. 

Hie  flux  of  CO2  from  the  soil  was  simply  estimated  using  measurements  of 
Moss  (1959)  as  a  guide. 

Successive  Approximations 

It  is  obvious  by  now  that  the  solution  to  each  system  of  equations 
presented  in  the  preceding  sections  of  this  chapter  has  depended  on 
the  solution  of  some  other  part  of  the  model.  This  interdependency 
makes  necessary  the  use  of  successive  approximations  to  solve  all  the 
equations  simultaneously.  At  the  beginning  of  this  routine,  the  soil,  air 
and  leaf  temperatures  were  set  equal  to  the  reference  temperature  at  the 
reference  height  above  the  crop.  The  air  water  vapor  pressures  were  also 
set  equal  to  the  value  at  the  reference  height.  Then  the  thermal  radi¬ 
ation  components  were  calculated.  Knowing  the  thermal  radiation  and  the 
air  temperatures  and  air  water  vapor  pressures  the  leaf  and  soil  temper¬ 
atures  were  calculated.  Then  the  sources  of  latent  and  sensible  heat 
were  calculated  which  were  used  in  turn  to  calculate  profiles  of  air  temp¬ 
erature  and  water  vapor  pressure.  Note  that  by  integrating  the  sources 
the  total  sensible  and  latent  heat  was  calculated  for  the  crop.  The 
procedure  was  repeated  until  the  values  of  the  sources  from  individual 
leaf  area  increments  converged  to  a  constant  value. 

A  similar  procedure  was  used  to  calculate  net  photosynthesis  once 
the  leaf  temperatures  were  known.  The  initial  values  of  the  C02  concen¬ 
trations  of  the  air  stream  were  set  equal  to  the  C02  concentration  at 
the  reference  height.  Knowing  the  air  C02  concentrations  sources  and 
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sinks  of  net  photosynthesis  were  calculated  which  were  in  turn 
♦-.J^.^te^calculate  a  CO^ profile.  .  .,  • .  •« 

Thus  the  complete  solution  resulted  in  not  only  the  source  and  sink 
distribution  of  net  photosynthesis  but  similar  distributions  of  sensible 
and  latent  heat  as  well.  Also  profiles  of  air  temperatures,  air  water 
vapor  pressures  and  CO2  concentrations  are  also  calculated.  These  prove 
extremely  useful  in  verifying  the  model  by  comparing  theory  with  experi¬ 


mental  results. 


Materials  and  Methods 

v  ,,  ^  «•  <  V  4  <•*  — -■**  ♦  *  ♦  •—  »  *  *’  •*  *• 

The  experimental  verification  of  a  mathematical  model  is  an  extremely 
important  part  of  its  development.  In  this  case,  field  measurements  of 
profiles  of  air  temperature,  COg  and  water  vapor  concentrations  along 
with  vertical  fluxes  of  CO^,  sensible  and  latent  heat  were  compared  to 
the  calculations  of  the  model.  Such  measurements  along  with  many  of  the 
necessary  plant  properties  (such  as  leaf  area,  leaf  angle  distribution, 
mesophjll  and  carboxylation  resistances,  etc.)  were  provided  in  a  com¬ 
prehensive  study  of  momentum  and  energy  exchange  above  and  within 
stands  of  corn.  The  complete  study  was  made  at  the  Microclimate  Experi¬ 
mental  Site  in  Ellis  Hollow  near  Ithaca,  New  York.  The  visible 'radia¬ 
tion  and  leaf  area  measurements  described  in  Chapter  I  were  part  of  this 
larger  study.  Hie  leaf  chamber  measurements  made  by  Dr.  R.  B.  Musgrave 
and  described  ->'n  Chapter  II  were  made  in  cooperation  with  this  study.  Hie 
complete  study  represented  contributions  from  a  team  of  individuals.  Hie 
author  was  the  team  member  directly  responsible  for  the  radiation  and  leaf 
area  measurements . 

The  field  measurements  were  made  in  a  20  acre  field  of  corn  which 
was  planted  at  a  density  of  6.0  plants/m^.  The  corn  (variety  -  Cornell 
M-3)  was  planted  in  38  cm  wide,  north-south  rows  with  approximately  46 
cm  between  plants  in  the  rows .  A  check  row  planter  was  used  to  place  the 
seeds  in  a  hexagonal  pattern.  Early  in  the  season  of  1968,  one-half 
the  field  was  thinned  to  4.4  plants/n£.  During  August  of  1968,  this 
same  portion  of  the  field  was  twice  more  thinned,  first  to  2.4  and  then 
to  1.4  plants/m^.  Measurements  in  each  half  of  the  field  were  made  before 


and.  after  each  thinning  in  August.  Measurements  in  the  thinned  and 

unthinned  portions  of  the  field  were  made  simultaneously. _  _ 

A  major  facet  of  the  field  study  consisted  in  applying  the  energy 
balance  to  estimate  vertical  fluxes  of  COg,  sensible  and  latent  heat. 
The  energy  balance  equation  can  be  written  as; 


RqCz)  =  LE(  z)  +  H(z)  +  -N  (z)  +  SHF 


(1UU) 


Net  radiation  (Rn)»  and  the  vertical  fluxes  of  sensible  heat  (H), 
latent  heat  (LE)  and  CC^N^)  are  expressed  as  functions  of  height  (z). 
Hie  constant,  X,  is  a  factor  which  converts  gms  of  CO^  into  energy 
units.  It  is  based  on  the  amount  of  energy  required  to  convert  a  mole 
of  CO2  into  a  mole  of  carbohydrate  (Lemon  (1962)).  The  above  fluxes  can 


be  expressed  as; 


LE(z)  =  -  ^(z)  3-eM 

Y*  9z 


(1^5) 


H(z)  =  -pCp  Kh(z) 


(1U6) 


Nt(z)  =  -X  Kh(z) 

0  n  9z 


(H»7) 


The  sign  convention  adopted  here  considers  upward  fluxes  as  positive. 
Thus  combining  equations  (1^5)  -  (l^7)  with  (lUU)  and  rearranging, 


results  in. 


K^z)  = 


Rn(z)  -  SHF 

_  9C(z)  _  c  T[zj_  _  PCP  3e(z) 
8z  P  8z  Y*  9z 


(1U8) 


Note  that  once  again  the  assumption  is  made  that  the  eddy  diffusivities 
for  mass  and  heat  transfer  are  identical.  Therefore  by  measuring  net 


radiation,  and  the  gradients  of  air  temperature,  water  vapor  and  C0o 
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concentrations  above  and  within  vegetation,  the  eddy  diffusivity  as  a 
function  of  height  can  be  calculated.  Then  by  substituting  the  values 
along  with  the  above  gradients  in  equations  (1U5)  -  (1U7)  the  vertical 
fluxes  of  COg,  sensible  and  latent  heat  can  be  calculated. 

Net  radiation  was  measured  in  the  stands  of  corn  with  the  same 
sampling  system  described  in  Chapter  I.  The  Fritchen  net  radiometer 
was  used  (Fritchen,  1965).  The  net  radiometer  simply  replaced  the  visible 
radiation  sensor  in  a  plexiglass  holder  and  was  pulled  back  and  forth 
through  the  crop  along  nine  meters  of  outrigging  cables  by  a  reversible 
motor  and  pulley  system.  Net  radiation  at  five  levels  in  each  stand  of 
corn  was  measured.  A  net  radiometer  was  also  mounted  above  each  stand 
of  com . 

Air  temperature,  and  C02  and  water  vapor  concentrations  were  measured 
using  a  complex  air  sampling  system.  This  system  was  designed,  con¬ 
structed  and  operated  by  M.  Johnson,  G.  Drake  and  Dr.  E.  R.  Lemon.  It 
will  be  described  here  only  in  very  general  terms.  Details  of  the  system 
can  be  found  in  Lemon  et  al.  (1970  )•  Bie  central  component  of  this  sampling 
system  was  a  double  shielded  air  intake  nozzle.  It  served  the  dual  role 
of  being  an  inlet  to  the  air  sampling  system  as  well  as  providing  an 
aspirated,  radiation  shielded  mount  for  thermocouples.  Ten  nozzles 
were  connected  in  parallel  at  a  particular  level  above  the  ground.  Air 
drawn  in  at  each  nozzle  was  sucked  into  a  manifold  which  was  in  turn  conn¬ 
ected  to  a  heated  copper  alloy  conduit  pipe.  This  pipe  transported  the 
air  from  a  given  level  to  laminated  polyethylene  -  aluminum  mylar  bags 
which  were  situated  in  an  air  conditioned  trailer  at  the  edge  of  the 
experimental  field. 


Each  nozzle  consisted  of  a  glass  tube  mounted  concentrically 
inside  a  larger  acrylic  tube.  The  thermocouples  were  concentrically 
mounted  inside  the  glass  tube.  For  radiation  shielding,  the  glass  tube 
was  covered  with  aluminized  mylar.  Also  the  acrylic  tube  was  covered 
with  a  1.25  cm  layer  of  styrofoam  which  was  in  turn  covered  with  alum¬ 
inized  nylar.  The  nozzle  had  an  outside  diameter  and  a  length  of  approxi¬ 
mately  5  and  25  cm  respectively.  One  thermocouple  in  a  nozzle  represented 
one  junction  of  a  ten  junction  thermopile.  Each  thermopile  was  referenced 
to  the  level  immediately  below  it.  Thus  there  were  two  sets  of  junctions 
at  each  level.  The  lowest  level  was  referenced  to  an  ice  bath. 

A  carbon-vane ,  oil- less  rotary  pump  drew  the  air  from  the  nozzles 
at  a  given  level  in  the  field  to  a  set  of  storage  bags  in  the  air  con¬ 
ditioned  trailer.  Air  was  sampled  simultaneously  at  all  levels  for 
one-half  hour  periods.  At  the  end  of  a  1/2  hour  period  the  air  from  all 
the  tubes  was  redirected  to  another  set  of  empty  bags.  Air  from  the  filled 
set  of  bags  was  then  sequentially  directed  through  an  infrared  water 
vapor  analyzer  (Modern  Controls)  which  measured  the  absolute  water 
vapor  concentration  at  each  level.  The  air  then  passed  through  the  sample 
cell  of  a  differential  infrared  COg  analyzer.  (Mine  Safety  Appliance 
Co.,  Model  LIRA  100 ).  Air  from  the  highest  level  above  the  ground 
was  directed  through  the  reference  cell.  The  CO2  concentration  at  this 
reference  height  was  also  measured  with  an  absolute  COg  analyzer.  (Mine 
Safety  Appliance  Co.,  Model  LIRA  200.)  The  millivolt  signals  from  the 
analyzers  were  recorded  on  strip-chart  recorders.  Converted  water 
vapor  and  CO2  concentrations  were  subsequently  plotted  versus  height 
above  the  ground  surface.  Gradients  were  estimated  from  smoothed  curves 
drawn  through  these  plotted  points. 


Wind  measurements  were  made  by  L.  H.  Allen,  Jr.  in  connection 
with  the  study  of  momentum  exchange.  Wind  velocity  was  measured  with  cup 
annemometers  above  the  vegetation ( Car di on  West,  C.  W.  Hiornthwaite 
Associates).  The  pulses  from  the  cup  annemometers  were  registered  on 
printing  counters.  (Machine  Electrification,  Model  MEK  ZDG  V).  These 
counters  printed  accumulated  counts  for  5  min.  intervals.  Within  the 
vegetation,  heated  thermocouple  annemometers  were  used  to  obtain  5  and 
30  minute  mean  values  (Hastings  -  Rpydist  Inc.,  Air  Meter  Model  Rm-IX, 
Probe  Model  N-7B). 

Incoming  total  and  diffuse  solar  radiation  and  total  and  diffuse 
visible  radiation  were  measured  as  a  part  of  the  study.  Solar  radiation 
was  measured  with  two  Epply  pyranometers .  Visible  radiation  was  measured 

with  two  custom  mounted  selenium  cells.  (The  design  of  this  sensor* is 

\ 

described  in  Chapter  i).  Shade  rings  were  placed  over  one  of  each  of 
the  above  types  of  instruments  to  record  the  incoming  diffuse  solar  and 
visible  radiation. 

Finally,  soil  heat  flux  plates  (National  Instrument  Laboratories, 
Inc.)  were  used  to  measure  the  heat  flux  into  the  soil.  A  set  of  plates 
were  placed  just  below  the  soil  surface  in  each  of  the  two  portions  of 
the  field. 

Millivolt  signals  from  the  heat  flux  plates,  all  of  the  radiometers, 
all  of  the  thermocouples  and  all  of  the  heated  thermocouple  annemometers 
were  fed  into  a  100  channel  data  logger  (A.  D.  Data  Systems  Incorporated). 
They  were  recorded  sequentially  on  a  7-track,  magnelic  tape.  (Digital 
Tape  Recorder:  Digi-Data  Corp.)  The  100  channels  could  be  sampled  as 
often  as  once  every  two  seconds.  The  7-track  tape  was  read  and  analyzed 
by  an  IBM  360/65  computer.  Average  half  hour  values  were  punched  out 


on  cards.  Values  from  the  smoothed  CO^  and  water  vapor  profiles  at 
small  height  intervals  were  also  punched  out  on  cards.  A  computer  program 
written  by  M.  Groom  read  the  values  of  air  temperature,  COg  and  water 
vapor  concentration,  net  radiation  and  the  soil  heat  flux  from  cards 
and  used  these  data  for  the  energy  balance  calculations.  The  program 
fit  the  air  temperature  and  net  radiation  data  to  polynomial  functions  of 
height.  The  temperature  polynomial,  when  differentiated,  supplied  the 
necessary  temperature  gradients.  Water  vapor  and  COg  gradients  were 
supplied  by  finite  differences  between  the  small  height  increments.  Thus 
the  program  was  able  to  calculate  the  eddy  diffusivity  as  a  function  of 
height  for  a  particular  set  of  data  along  with  the  vertical  fluxes  of 
CC>2 ,  sensible  and  latent  heat. 

A  study  of  plant  water  relations  was  made  simultaneously  with  the 
above  measurements  by  R.  W.  Shawcroft.  Leaves  near  the  top  of  the  crop 
and  near  the  bottom  were  sampled  periodically.  The  leaf  stomatal  re¬ 
sistance  using  a  Van  Bavel  porometer  (Van  Bavel  e£  al .  1965 )  and  leaf 
relative  water  content  were  measured.  The  soil  moisture  tension  using 
tensiometers  was  also  measured  in  connection  with  this  water  relations 
study. 

The  above  measurements  represented  a  set  of  parameters  which 
provided  information  needed  in  the  model  as  well  as  numbers  with  which  to 
verify  the  model.  At  this  point  a  summary  of  information  required  by 
the  model  is  in  order. 

First  of  all  the  values  of  air  temperature,  CO2  and  water  vapor 
concentration,  and  wind  velocity  are  needed  at  some  reference  height  above 
the  crop.  Also  values  of  net  radiation,  direct  and  diffuse  solar  radiation 


and  direct  and  diffuse  visible  radiation  are  needed  above  the  crop. 

Finally  the  soil  heat  flux  and  soil  moisture  tension  at  the  soil  surface 
as  well  as  the  soil  albedo  are  needed. 

However,  the  above  mentioned  boundary  conditions  are  only  part  of 
the  input.  The  model  itself  is  a  complex  function  of  plant  characteristics 
which  have  to  be  known  as  well  before  any  solution  can  be  obtained.  These 
plant  parameters  are  of  two  types. 

The  displacement  plane  D,  the  roughness  length  zQ  and  the  adjusted 
crop  height  (CAHT)  which  have  been  discussed  earlier  in  this  chapter  are 
properties  of  the  plant  community  as  a  whole. 

On  the  other  hand,  individual  leaf  characteristics  such  as  leaf 
transmittance  and  reflectance  as  described  in  Chapter  I  as  well  as 
mesophyll  and  carboxylation  resistance  and  stomatal  resistance  as  a 
function  of  light  intensity  are  also  needed. 

Finally  the  leaf  area  density,  the  cumulative  leaf  area  index  and 
the  leaf  angle  distributions  as  described  in  Chapter  I  represent  a  most 
important  input  to  the  model  because  we  are  particularly  interested  in 
studying  the  effects  of  crop  architecture  on  the  various  exchange  processes 
already  mentioned.  Only  by  incorporating  a  detailed  description  of  the 
plant  community  in  the  mathematical  model,  can  one  do  this. 
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Results  and  Discussions 

At  this  time,  the  crop  simulation  model  has  not  been  extensively 
tested  against  the  measurements .  Three  test  periods  on  the  1968  unthinned 
com  (1145-1215,  Aug.  15;  0845-0915,  Aug.  18  and  1145-1215,  Aug.  18)  are 
reported  here  as  examples .  The  test  periods  were  relatively  cloud  free 
so  that  incoming  solar  radiation  was  almost  constant  over  each  half-hour 
interval.  However,  the  test  periods  differed  in  that  August  15  was 
much  drier  than  August  18.  The  drier  conditions  were  clearly  reflected 
in  the  leaf  resistance  and  water  content.  The  effect  of  water  stress 
will  be  considered  in  detail  below. 

First  of  all  it  is  important  that  wind  velocity  can  be  predicted  in 
the  vegetation  from  measurements  at  the  reference  height.  The  wind 
velocity  is  used' in  equation  (124)  to  calculate  boundary  layer  resistances. 
Figure  (30)  represents  a  typical  wind  profile  with  D  =  140  cm,  z0  =  IT  cm 
and  CAHT  =  180  cm.  This  profile  was  measured  in  the  period  1145-1215, 
August  15.  In  general,  for  the  purposes  of  the  model,  the  theory  pre¬ 
dict..  *  the  wind  velocity  with  height  in  the  crop  quite  adequately.  There 
is  some  discrepency  at  the  bottom  of  the  crop.  Here  ■‘■he  measured  values 
remain  constant  while  the  theoretical  values  continue  to  decrease.  This 
discrepancy  is  not  serious  to  the  model  output  but  is  interesting  nonethe¬ 
less.  The  measurements  indicate  an  entrainment  of  air  underneath  the 
vegetation. 

The  eddy  diffusivities  for  heat,  calculated  by  the  model  show  a 
number  of  differences  from  the  values  calculated  from  energy  balance 
measurements  as  illustrated  in  Figure  (  3l)  •  First,  the  model  consistently 
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Figure  3l(a).  Eddy  diffusivity  for  sensible  heat  (Kh)  as  a  function  of 
height  (s)  (Aug.  15,  1200). 


Figure  3l(b).  Eddy  diffusivity  for  sensible  heat  (K  )  as  a  function  of 
height  (z)  (Aug.  18,  1200).  1 
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Figure  bl(c). 


Eddy  diffusivity  for  sensible  heat  (%)  as  a  function  of 
height  (z)  (Aug.  18,  0900). 


underestimated  the  eddy  diffusivities  above  the  crop.  In  the  crop,  the 
theoretical  values  of  decrease  approximately  in  an  exponential 
manner.  However,  the  measured  mid-day  profiles  of  Kjj  in  the  crop  show 
a  curious  "S"  shape  as  one  moves  down  into  the  crop.  This  'dip'  in  the 
profile  is  not  as  apparent  as  the  0845-0915  test  period.  The  dip  remains 
a  mysw  y  at  this  time.  Further  examination  of  more  data  is  in  order 
before  making  any  suggestions  as  to  its  cause. 

Figures  (32-34)  shows  theoretical  and  measured  profiles  of  CO2  and 
water  vapor  concentrations,  air  temperature  and  net  radiation.  The 
total  vertical  fluxes  above  the  crop  of  CC^  and  sensible  and  latent  heat 
are  given  in  Table  II.  Although  the  agreement  between  theory  and  meas¬ 
urements  if  far  from  perfect  there  are  consistent  trends  in  the  data. 

The  most  obvious  difference  is  the  smaller  amount  of  sensible  heat  and 
the  larger  amount  of  latent  heat  generated  by  the  model. 

This  larger  amount  of  sensible  heat  could  explain  in  part  the  larger 
eddy  diffusivities  above  the  crop.  The  larger  sensible  heat  flux  increases 
the  z/L  ratio  which  was  described  in  the  theory.  This  in  turn  increases 
the  Kjj/K^  ratio.  However,  even  on  August  15  when  the  theory  and 
measurements  came  closest  to  agreeing,  measured  by  the  energy  balance 
method  was  approximately  twice  as  large  as  the  theory  predicted.  An 
unsteady  state  condition  in  the  field  could  also  lead  to  a  larger 
exchange  coefficient.  This  larger  eddy  diffusivity  from  the  measurements 
made  matching  theoretica.1  to  measured  profiles  much  more  difficult.  Thus, 
boundary  values  in  the  model  were  set  at  values  which  matched  the  theor¬ 
etical  values  of  air  temperature,  COg  and  water  vapor  concentration  to 
measured  values  at  the  top  of  the  crop.  This  made  comparison  within  the 
crop  much  easier  in  Figures  (32-34). 
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Figure  33(a).  Theoretical  profiles  of  C02  and  water  vapor  concentration 
(e)  with  measurements  (Aug.  18,  1200).  The  CO^  concen¬ 
trations  as  shown  are  subtracted  from  the  reference  height 
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Figure  3U ( a. ) .  Theoretical  profiles  of  C<  o  and  -rater  vapor  concentration  (e) 
with  measurements  (Aug.  16,  090j).  The  C02  concentrations 
as  shown  are  subtracted  from  +he  concentration  at  the 
reference  height. 
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Mean  theoretical  and  energj  oala.  ::e  'lues  for  the  total  flux  of  sensible 
and  latent  heat  and  of  COg  J’Vo~  tne  unthinned  corn  crop. 
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0.51+ 
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0.71 


.025 
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Immediately  apparent  in  the  within  the  crop  comparison  is  the 
relatively  large  difference  between  the  theoretical  and  measured  air 
temperatures.  Although  the  theoretical  water  vapor  concentrations  are 
in  general  higher  than  the  measured  in  keeping  with  the  higher  theoreti¬ 
cal  amounts  of  latent  heat,  the  discrepancy  is  not  as  apparent  as  in  the 
temperature  profile. 

These  temperature  differences  in  the  axi  stream  could  be  caused  by 
an  inadequate  theoretical  treatment  at  the  soil  surface.  There  is  a 
stronger  source  of  sensible  heat  in  the  field  than  predicted  by  the  model. 
This  extra  heat  could  come  in  part  from  the  soil  surface  itself.  Ihe 
visible  radiation  study  of  Chapter  I  indicated  that  approximately  7-10JS 
more  radiation  penetrates  to  the  soil  surface  than  is  predicted  by  the 
model.  The  soil  surface  itoelf  at  the  experimental  site  is  to  a  large 
extent  covered  with  stones  which  are  dry  and  heat  up  readily  since  rock 
has  a  lower  specific  heat  than  moist  soil.  Also  the  latent  heat  flux  at 
the  soil  surface  is  very  difficult  to  calculate.  This  arises  from  the 
difficulty  in  predicting  mean  or  effective  soil  moisture  tension  at  the 
surface.  In  the  model,  a  value  is  assumed  which  gives  approximately  equal 
sensible  and  latent  heat  fluxes  from  the  soil  surface.  This  combination 
of  factors  leads  to  an  underestimation  of  the  soil  surface  temperature. 

It  is  significant  that  the  discrepancy  between  theory  and  measurements  is 
less  at  0900,  Aug.  18  when  the  amounts  of  radiation  reaching  the  soil 
surface  is  about  one-half  of  the  amount  at  1200  hours. 

The  model  does  show  good  agreement  between  theory  and  measurements 
both  in  the  profiles  of  CO2  and  in  the  estimations  of  net  photosyr thesis. 
There  is  a  small  consistant  overestimation  by  the  model.  However, 


considering  the  errors  involved  in  the  energy  budget  measurements  and  the 
uncertainties  associated  with  stomata!  resistance  and  leaf  chamber 
measurements  along  with  the  various  assumptions  made  in  the  model,  the 
agreement  between  theory  anc’.  measurement  for  net  photosynthesis  is  quite 
adequate . 

What  is  especially  noteworthy  in  these  test  periods  are  the  resulting 
effects  of  increasing  stomatal  resistance.  The  averaged  measured  leaf 
resistance  in  the  uppermost  leaves  at  1200  on  August  15  was  2.87  sec/cm. 

This  average  value  was  very  close  to  0.97  for  the  time  periods  of  August  18. 
Thus  for  the  dry  day,  y0,  the  minimum  leaf  resistance  was  changed  from 
0.97  to  2.87  sec/cm.  This  was  reflected  in  the  theoretical  amounts  of 
sensible  and  latent  heats  which  were  much  closer  to  the  measured  values  for 
this  test  period.  Curiously,  the  energy  balance  measurements  of  latent 
heat  and  sensible  heat  do  not  reflect  the  effects  of  larger  stomatal  re¬ 
sistances  on  the  drier  day  (Aug.  15). 

The  energy  balance  measurements  do  agree  with  the  model  when  estimat¬ 
ing  net  photosynthesis  in  that  both  theoretical  and  measured  values  of 
net  photosynthesis  are  higher  on  the  day  of  smaller  leaf  resistances 
for  the  same  time  periods  (1200).  At  this  time,  both  days  had  approxi¬ 
mately  equal  amounts  of  solar  radiation.  Air  temperature,  and  wind 
velocities  were  similar  as  well.  While  there  are  many  other  factors 
involved,  there  is  evidence  here  to  suggest  that  the  leaf  resistance  had 
a  significant  effect  on  the  photosynthesis  rate.  However,  more  extensive 
comparisons  are  needed  to  verify  the  effect. 
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Conclusions 

The  limited  comparison  of  the  model  to  field  measurements  made  in 
this  study  is  not  meant  to  be  a  critical  test  of  the  model.  Much  more 
extensive  work  remains  before  the  model  can  be  truly  evaluated.  How¬ 
ever  ,  in  this  initial  stage,  the  model  does  show  signs  of  becoming  a 
powerful  tool  to  be  used  in  the  understanding  of  crop-environment  inter¬ 
actions.  Especially  encouraging  has  been  the  manipulation  of  the  stomatal 
resistance  to  give  a  decrease  in  photosynthesis.  This  provides  a 
means  to  study  parameter  significance  in  the  water  use  efficiency  of 
plants . 

Improvements  over  previous  models  are  that  no  extinction  coefficient 
had  to  be  assumed  for  the  penetration  of  visible  or  net  radiation  into 
the  crop  as  long  as  the  distribution  of  the  leaf  area  is  random.  Also 
no  extinction  coefficient  was  used  to  describe  the  eddy  diffusivity  with 
height  in  the  crop.  The  eddy  diffusivity  was  coupled  directly  with  the 
wind  velocity  which  was  in  turn  predicted  in  the  crop.  The  treatment 
of  visible  radiation  distribution  on  plant  leaves  was  extended  to  include 
the  distribution  of  near  infrared  radiation  as  well.  This  along  with  a 
treatment  of  thermal  radiation  allowt  1  the  net  radiation  to  be  calculated. 

This  model  also  represents  an  attempt  to  incorporate  the  effects 
of  leaf  area  density  on  the  mixing  processes  in  the  crop.  However,  in 
general  the  theory  did  not  adequately  reflect  the  energy  balance  measure¬ 
ments  both  above  and  within  the  crop.  A  much  fuller  treatment  of  turbulent 
exchange  in  vegetation  still  needs  to  be  developed. 

Similarly  there  is  still  much  work  to  be  done  in  effectively  dealing 
with  the  soil  surface.  It  is  especially  difficult  to  estimate  the  flux  of 


water  from  any  soil  surface  because  of  the  difficulty  in  estimating  the 
soil  moisture  tension.  More  detailed  models  will  include  treatments  of 
unsaturated  water  flow  in  soils  to  the  soil  surface  and  to  roots.  The 
model  by  Cowan  (1965)  is  a11  example  of  this  type  of  work. 

In  general,  the  use  of  mathematical  models  helps  considerably  in 
systematically  attacking  the  problems  of  understanding  the  complex 
interactions  between  the  plant  community  and  the  environment.  This  model 
attempted  to  sort  out  and  evaluate  factors  which  contribute  to  net 
photosynthesis  over  a  short  time  period.  While  this  in  no  way  represents 
an  understanding  of  all  the  growth  processes,  it  does  allow  one  aspect 
of  growth  to  be  singled  out  and  studied  intensively. 
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